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ABSTRACT 

Techniques  for  extrapolating  statistical  data  on  the  parameters 
of  individual  components    to  the  statistical  performance  index  for  an 
overall  system  are  considered.      Two  cases  are  evaluated.      (1)    The 
deterministic  case  in  which  the  system's  performance  index  is  known 
functionally  in  terms  of  the  system  parameters.      (2)     The  non- 
deterministic  case  in  which  only  limited  data  on  the  performance  index 
and  its  sensitivity  with  respect  to  system  parameters  is  known.      Com- 
puter programs  are  developed  in  both  cases  for  combining  given 
probability  density  distributions  for  the  parameters  into  an  overall 
probability  density  distribution  for  the  system  performance  index. 
Theory  and  programs  are  developed  and  verified  with  a  specific  numer- 
ical example. 
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I.       INTRODUCTION 

The  performance  of  a  system,    from  a  statistical  point  of  view, 
depends  upon  the  nominal  values  of  the  parameters  of  the  component 
parts  included  in  the   system  and  on  the  statistical  distributions  of 
those  parameters.     Alternately,    the   response  of  the  component  parts 
to  a  given  signal  may  be  known  statistically  and  it  is  desired  to  pre- 
dict an  overall  statistical  distribution  for  the   system  response. 
Statistical  information  on  the  system's  response  may  be  interpreted 
in  terms  of  probability  of  success  or  failure  (operation  or  non- 
operation)  for  a  single  system  or  in  terms  of  the  percentage  of 
operational  systems  that  can  be  expected  in  an  ensemble  of  systems. 
Currently  three  basic  methods  are  being  used  to  obtain  performance 
data  for  systems;  the  worst-case  method  [Mark  and  Stember  1965], 
the  moment  method  [Mark  and  Stember  1965],    and  the  Monte  Carlo 
method  [Shreider  1966]. 

The  worst- case  method  is  a  non- statistical  approach  intended  to 
determine  whether  it  is  possible,    with  given  parameter  tolerance 
limits,    for  the  systems  performance  characteristics  to  fall  outside 
specifications.     This  method  determines  nothing  of  the  probability 
distribution  of  the  performance.      Worst-case  analysis,    although  it 
may  be  involved,    does  provide  useful  engineering  insight  and  is  used 
a  great  deal  in  practice.      However,    its  outlook  is  basically  limited. 


Statistics  are  combined  with  system  analysis  techniques  in  the 
moment  method  to  estimate  the  probability  that  performance  will  re- 
main within  specified  limits.      The  actual  parameter  probability 
distributions  are  usually  replaced  by  Gaussian  distributions  having 
the   same  variances.      Hence,    the  moment  method  has  inherent  limita- 
tions.     Parameter  probability  distribution  characteristics  such  as 
skewness  and  peakedness  are  usually  not  taken  into  account. 

In  the  Monte  Carlo  method  a  large  number  of  alternate   replicas 
of  a  system  are  simulated  by  mathematical  models.      Component 
values  are   randomly  selected  and  the  performance  of  each  replica  is 
determined  for  its  particular  set  of  components.      The  accuracy  of  the 
Monte   Carlo  method  is  dependent  upon  taking  a  sufficiently  large 
number  of  replicas  and  can  be  both  lengthy  and  costly.      Monte   Carlo 
methods  are  also  used  extensively  in  practice  although  they  provide 
limited  engineering  information. 

The  purpose  of  this  thesis  is  to  present  a  numerical  technique 
for  combining  statistical  data  on  several  independent  components  into 
a  single  distribution  function  for  the   system  performance.      This  is 
accomplished  by  considering  two  general  cases.      The  deterministic 
case,    where  the  functional  relationships  between  the  system  parameters 
and  the  output  are  known  functionally;  and  the  non-deterministic  case, 
where  these  functional  relationships  are  unknown  but  measured  dis- 
crete data  is  available.      In  both  cases  the  actual  parameter  prob- 
ability distributions  are  used  in  conjunction  with  numerical  integration 
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techniques  to  obtain  the  system  distribution  function.      One  major 
advantage  of  the  approach  is  that  the  known  parameter  distributions 
are  quite  arbitrary  and  may  be  histograms  obtained  from  the  results 
of  experiments  on  individual  components. 
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II.       DETERMINISTIC  CASE 

The  object  of  the  deterministic  case  is  to  investigate  the  dis- 
tribution of  a  system  output  when  the  output  has  a  known  mathematical 
relationship  to  several  internal  parameters.      The  basic  assumptions 
are  that  all  of  the  internal  parameters  are  independent  and  that  the 
probability  distributions  of  the  internal  parameters,    although  arbi- 
trary,   are  known  mathematically  or  given  in  terms  of  histograms. 

To  commence  the  investigation,    the  general  equations  for  various 
combinations  of  parameter  distributions  which  may  be  expected  in  a 
typical  system  performance  index  must  be  obtained.      These  will 
basically  consist  of  operations  on  the  distributions  of  random  variables 
which  are  interrelated  by  addition,    subtraction,    multiplication,    and 
division.      The  basic  mathematics  for  these  operations  are  well  known 
in  the  literature,    at  least  in  general  form  [Beckmann  1967], 

Consider  two  independent  random  variables:     X,    assuming  values 
x  with  a  known  probability  density  py^(x),    and  Y,    assuming  values  y 
with  a  known  probability  density  py(y).     Also  consider  the  functions 
U  =   U(X,  Y)  and  V  =  V(X,  Y)  and  their  inverse  functions  X  =  X(U.  V) 
and  Y  =  Y(U,  V).      Then,    when  X  is  near  x  and  Y  is  near  y,    U  and  V 
must  be  near  u  =  u(x,  y)  and  v  =  v(x,  y).      Hence, 

PXY(x,  y)  dxdy  -  puv(u,  v)  dudv  (1) 

Substituting  the  inverse  functions  and  the  required  probability 

densities,    we  obtain 
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PTTV(u,  v)  =  p       [x(u,  v)]  pv  [y(u,  v)] 


uv 


'X 


Y 


i(u,v) 


(2) 


where  pxjy(u,  v)  is  the  joint  probability  density  of  U  and  V,    p^-[x(u,  v)] 
is  the  probability  density  of  X,    py  [y(u,  v)]  is  the  probability  density 


of  Y,    x(u,  v)  and  y(u,  v)  are  the  inverse  functions,    and 
is  the  Jacobian,    given  by 


l(u,v) 


3X 


ax 
iZ 


(3) 


When  X  and  Y  are   related  by  a  single  functional  the  foregoing  can  be 
modified  by  letting     U  =  U(X,  Y)  and  V  =  X.  (4) 

Then, 


dX 

1 

au 

„  ,.  HL 

*u 

M 

0 

du 

a(x<,.y 
d(u,v 

Equation  (2)  then  becomes 

Pux(u,  x)   =  px(x)py  [y(u,  x)] 


(5) 


iZ 


au 


(6) 


r00 

Since,    in  general,    prj(u)     =        prr-vr  (u,  x)dx,    equation  (6)  may  be  written 


as 


Pu(u)  =    I  px(x)py  [y(u'x>]  i^l dx 

-oO 


(7) 


Using  equation  (7),    the   specific  formulas  of  Table  I  have  been 
compiled  for  the  functional  relationships  and  operations  listed,    which 
were  chosen  on  the  basis  of  what  might  be  expected  in  engineering 
problems. 
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The  technique  outlined  here  consists  of  the   successive  application 
of  the  formulas  of  Table  I  as  required,    with  the   required  integration 
performed  numerically.      Various  types  of  numerical  integration  are 
possible;   rectangular  integration,    trapezoidal  integration,    or  Simpson's 
Rule.      These  are  summarized  as  follows- 

Rectangular  integration  [Stanton  1961], 


Xn  *1 


Trapezoidal  integration  [Stiefel  1963], 


(8) 


/Xin 


(9) 


(x)  AxJ2[f(xl)+f(xn)]  45Zf(xi)l 

'*•  ^  i=2 

Simpson's  Rule  [Stanton  1961], 
,-X6+nAX 

J  A  YL  (10) 

where     C    =     1,    4,    2, 2,    4,    1 

Although  the  specific  technique  of  numerical  integration  is  of 
prime  importance  in  so  far  as  accuracy  and  computer  programing 
and  time  are  concerned,    the  exact  procedure  used  is  a  detail  which 
is  best  dictated  by  the  particular  problem  to  which  these  techniques 
are  to  be  applied. 

This  method  may  be  best  illustrated  by  considering  the  following 
example.      Consider  that: 

H    =  (11) 
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where  a  and  b  are  constants  and  X,  Y,    and  Z  are  independent  random 
variables  with  known  (measured)  probability  densities.      Consider  that 
p-^(x),    pY(y),    and  p7(z)  are  given  as  follows: 


(a)     Consider  first  the  combination  U  =  XY,    with  4X     =     0.  1. 
Using  rectangular  integration  for  simplicity,    equation  (3)  of  Table  I 
becomes,    for     u  =  U-     : 


^i>  =j^<*)PY<»/s>tr  =  **f-*~ 


PliiT 


When     u  =  u. 

=  0.  09 

-Xi 

0.9 

0.  1 

4.  0 

0.  2 

0.  45 

4.  0 

0.  3 

0.  3 

4.  0 

0.  4 

0,  22 

4.  0 

0.  5 

0.  18 

4.  0 

0.  6 

0.  15 

2.  6 

0.  7 

0.  13 

1.  7 

0.  8 

0.  11 

1.  0 

0.  9 

0.  1 

0.  44 

i=l 


APfa)     ¥ 

"M«M 

1.  8 

7.  2 

0.9 

3.  6 

0.6 

2.  4 

0.  43 

1,  72 

0.  37 

1.  48 

0.  30 

0.  78 

0.  26 

0.  44 

0.  21 

0.  21 

0.  20 

0.  088 

sum  = 

17. 918 

12) 


Putu^  =  (0.  1)(17.  918)  =  1.79 
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when   u=u?= 

0.  2 

Xi 

U*A; 

0.  1 

2.  0 

0.  2 

1.  0 

0.  3 

0.  7 

0.  4 

0.  5 

0.  5 

0.  4 

0.  6 

0.  33 

0.  7 

0.  3 

0.  8 

0.  25 

0.  9 

0.  22 

0.  0 

-V*A 

4.  0 

0.  0 

4.  0 

2.  0 

8.  0 

4.  0 

1.4 

5.  6 

4.  0 

1.  0 

4.  0 

4.  0 

0.  8 

3.  2 

2.6 

0.  67 

1.  74 

1.  7 

0.6 

1.  02 

1.  0 

0.  5 

0.  5 

0.44 

0.  42 

0.  185 

sum    =      24.  245 
Pu(u2)  =  (0.  1)  (24.245)  -  2.4 

The  pattern  for  evaluating  Ptt(u)  is  repeated  for  u  =  .  3,    .4,    --,    1.0 
yielding  the  curve  of  Figure  1. 


(b)     Next  consider  the  function  V 
(7)  of  Table  I  may  be  applied. 


1_ 

U 


XY 


Here  equation 


(13) 


— 
0.1 

10.  0 

1.  8 

_ 
0.  01 

0.  018 

0.  2 

5.  0 

2.42 

0.  04 

0.  097 

0.  3 

3.  33 

2.  39 

0.  09 

0.  215 

0.  4 

2.  5 

2.  11 

0.  16 

0.  338 

0.  5 

2.  0 

1.  54 

0.  25 

0.  385 

0.6 

1.  66 

1.  02 

0.  36 

0.  364 

0.  7 

1.43 

0.  55 

0.49 

0.  270 

0.  8 

1.  25 

0.  28 

0.  64 

0.  179 

0.9 

1.  11 

0.  05 

0.  81 

0.  040 

Py(v)    is  plotted  in  Figure  2. 
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(c)    The  function    W    may  now  be  calculated  using  equation  (8) 

of  Table  I,    where  W    =     aV  +  b=^_+b  (14) 

.X.  Y 


and 


-ri      °  |al £  v 


|alpVw 


(15) 


If    a    =     2  and  b    =     1,    the  following  result  is  obtained. 


V 

W  =  aV  +  b 

0.  0 

pv(v)/2 

1 

3 

0.  0 

2 

5 

0.  385 

0. 1925 

3 

7 

0.  25 

0.  125 

4 

9 

0.  15 

0.  075 

5 

11 

0.  10 

0.  050 

6 

13 

0.  06 

0.  030 

7 

15 

0.  04 

0.  020 

8 

17 

0.  03 

0.  015 

9 

19 

0.  02 

0.  010 

10 

21 

0.  01 

0.  005 

The  result  is  plotted  in  Figure  3. 


(d)    Since,    from  equations  (11)  and  (14), 


H    =     W    +     Z 


(16) 


it  follows  from  equation  ( 1 )  of  Table  I  that  Ppj(h)  is  given  by  the  con- 
volution of  pJw)  and  p„(z).      This  operation  may  also  be  carried  out 
numerically  and  for  the  given  p7(z)  may  be  expressed  as 


Pjj(h)  I  for 


PH 


(V.)  fo: 

h-1 


£  hs=l 


h  >  1 


17) 


These  integrals  may  also  be  carried  out  numerically  yielding  the 
results  plotted  in  Figure  4, 
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To  work  more  complicated  problems  than  the  previous  example 
requires  a  great  amount  of  mathematical  calculations.      Thus,    when 
using  numerical  integration  techniques  the  method  lends  itself  nicely 
to  the  digital  computer.      Figures   5  through  9  are  the  computer  flow 
graphs  corresponding  to  some  of  the  equations  of  Table  I.      The  com- 
puter program  for  any  mathematical  function  can  then  be  derived 
from  the  appropriate  combinations  of  these  basic  flow  graphs. 

To  illustrate  a  more  complex  problem  consider  the  following 
example- 


given 


i) 


j-       - 


where 


A    =     B3aE6a  .         B    =     B3bB6b 


and 


B.  B3b 


3a 


1     +     B3b(2.  16  x  10-2) 


B         =  B6b 

6a  1     +     B6b(0.63  x  10-2) 

B_,      is  a  normal  distribution;    mean  =  130,    standard  deviation  =  10 
3b 

B^b    is  a  normal  distribution;    mean  =  65,    standard  deviation  =   12 
W  is  a  uniform  distribution  from  w  =  -  4/7  to     w  =  +  4/7 

Theta  is  system  output  parameter  and  the  variables  B^,    B/v  ,    and 
W  are  independent. 
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Since  the  variables  are  independent,    equation  (18)  can  first  be 
simplified  by  ordinary  mathematical  manipulation. 

Let    X     =     B  Y    =     B6b,    kj_     =     2.  16  x  10"2,    k2  =  0.  63  x  10"2 

then  B    =    XY 


A     = 


XY 


(l     +    ktX)  (1     +    k2Y) 


1 
A 


1 
B 


and    8    =     0.  175    +      27'  8        + 

Y 


klk2    + 


8.  1 


X 


X 


w 


(i9: 


Now,    utilizing  the  basic  methods  and  flow  graphs,    the  problem 
can  be  programed  on  the  digital  computer  and  the  probability  density 
of  9     readily  obtained.     The  program  and  solution  are  contained  in 
Appendix  I. 
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III.        NON-DETERMINISTIC  CASE 

In  many  practical  problems  the  functional  relationship  between 
variables  is  unknown  and  it  cannot  be  expressed  simply  by  formula. 
It  is  generally  available  only  as  numerical  data  from  a  computer  run 
on  a  mathematical  model  of  the  system  under  study.      Two  cases  are 
considered  here.      First,    the  case  where  the   system  response  of 
interest,    U(X),    and  its  sensitivity  coefficients,    U'(X),    are  available 
at  discrete  values  of  a  single  variable,    X. 


U(X)       =      [u(Xl),    u(x2), ,    u(xn)] 

(X)     =      [u'fx^,    u'(x2), ,    u'(xn)] 


U 

Second,    the  case  where  the   system  response  of  interest,    U(X,  Y), 
and  its  sensitivity  coefficients,    U'(X,  Y),    are  available  at  discrete 
values  of  two  independent  variables,    X  and  Y. 


(20) 
(21) 


U(X.  Y)     = 


Uii      U,  -.       -    -    -    u 


where 


and 


11     u12 
u21     u22 

lml  um2 


uij     =     u(x.,  y.) 


In 

<sn 


mn 


(22) 
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U'(X,  Y) 


u'll      u'l2 


u  21      u  22 


ml         m2 


u 


-  u' 


In 


2n 


u 


t 
mn 


(23) 


where 


u' 


ij 


J>u(xi  ,    y.) 

ex- 


In  both  cases  the  problem  reduces  to  determining  Ptt(u)  from  the 
known  distributions  py(x),    Pv(y)  and  the  data  of  equations  (20)  and  (21), 
or  (22)  and  (23). 

In  the  first  case,    consider  that  the  probability  density,    p,,(x), 
and  the  cumulative  distribution, 


(24) 


are  known.      Then  at  each  data  point  of  u(x), 


U 


[  u(xj»] 


px  <xj> 


(25) 


That  is,    the  probability  that    u    ^    u(x.)  is  equal  the  probability  that 
x   —    x.  .      Equation  (2  5)  gives  discrete  points  on  the  cumulative 
distribution,    Ptt(u),    for  each  value  of  the  parameter  x  for  which 
u(x)  has  been  measured. 

It  is  also  well  known  [Beckman  1967]  that 


Pu(u) 


A  U 


' 


(26) 
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I 

- 

Thus  the  slope  of  the  cumulative  distribution  curve,    - = 

at  each  data  point  is  equal  to  the  known  probability  density,    py-(x-), 

divided  by  the  measured  sensitivity  function,  ,    at  the 

data  point.      The  conclusion  is  that  the  cumulative  distribution,    PTT(u), 

and  its  slope  may  be  calculated  readily  from  the  measured  response 

and  sensitivity  data.      This   curve  is  therefore  evaluated  as  a  piece - 

wise  linear  approximation. 

When  U  is  a  function  of  two  variables,    X  and  Y,    the  foregoing 

cumulative  distribution  curve  can  be  evaluated  for  each  value  of  the 

second  variable,    Y.      Thus  a  family  of  such  curves  or  data  points 

denoted  by  Ptj(u,  y-)  is  available.      For  a  specific  value  of  u  =  u- 

each  of  these  curves  gives  the  probability  that  u  ^  u-    when  y  =  y-. 

The  total  probability  that  u -^  u-  depends  upon  Py^AV       and  is 

given  by  the  equation: 

m 
P^Ui)     =    Ayy"Pu(ui>yi  j)  (27) 

.1=1 

where  Ptt(u,  y)  is  the  cumulative  distribution  evaluated  with  y  as  a 
constant  and  x  as  a  single  variable  as  discussed  previously. 

This  result  can  be  proven  by  considering  the  following  (see 

equation  (7)  ): 

f  °°  -t 

(28) 


Pu(u)   =    K('yjpx[x'(u,y)]|> 

-00 


where  u  =  u(x,  y).  is  the  sensitivity  function  with  respect  to 

the  parameter  x  and  in  equation  (28)  must  be  expressed  in  terms  of 
u  and  y.      It  is  convenient  to  denote 
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(29) 


and  then  equation  (28)  becomes 


(30) 


We  also  know 


U)        =  PyCu) 


(31) 


Substituting  equation  (30)  into  equation  (31)  and  inter- changing  the 
order  of  integration  yields 

J  J  .  A 


Pu(u) 


(32) 


Now 


itfu  r,  PTFY(u,y) 

-oo  -oo 

Since,    if  Y  is  considered  a  constant  and  X  as  the  random  variable, 


u  =  u(x,  y), 


pu<u'  y) 


px     [x(x,  y)]  px     [x(x,  y)] 


(34) 


AH 
dx 


ux(u,  y) 


Substituting  equation  (33)  into  equation  (32)  yields: 


-    OO 


Pn(u)     =         I   pY(y)Pu(u,  y)dy 


(35) 


In  discrete  form,    equation  (35)  becomes  equation  (27)  when 
rectangular  integration  is  used.      However,    more  accurate  integration 
techniques  may  also  be  applied  to  equation  (3  5). 
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To  demonstrate  the  application  of  this  non-deterministic  method 
consider  part  of  the  previous  example  of  equation  (18).      It  is  now 
assumed  that  the  functional  relationship  between  the  output  and  the 
parameters  is  unknown.      However,    we  do  know  the  probability  den- 
sities of  the  parameters  x  and  y  and  have  certain  measured  data  which 
consists  of  measured  values  of  x,    y,    and  0  and  the  corresponding 
sensitivities  for  each  measured  0.     The  data  is  presented  in 
Table  II.      Following  the  procedure  outlined,    the  problem  can  be 
programed  on  the  digital  computer  and  the  probability  density  of  0 
readily  obtained.     The  program  and  solution  are  contained  in 
Appendix  II.     The   results  can  be  compared  with  those  obtained  using 
the  deterministic  formulas.      Noting,    of  course,    that  the   solution 
to  the  deterministic  problem, contained  in  the  Appendix,    includes 
the  addition  of  the  uniform  distribution  W(w)  and  the  solution  to 
the  non-deterministic  problem  does  not. 
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IV.    CONCLUSIONS 

The  numerical  technique  presented  here  has  distinct  advantages 
over  the  worst-case,    the  moment,    and  the  Monte  Carlo  methods.      By 
utilizing  the  actual  parameter  probability  distributions  the  results 
obtained  are  more  accurate  than  those  of  the  moment  method  and  pre- 
sent a  great  deal  more  information  than  the  worst- case  method.      The 
accuracy  obtained  is  equivalent  to  that  of  a  large  number  of  Monte 
Carlo  replicas.     Thus  the  saving  in  time  over  the  Monte  Carlo  method 
is  obvious.      In  general,    the  numerical  method  presented  here  can 
obtain  a  more  accurate  solution  in  less  time  than  any  of  the  other 
current  methods.      However,    certain  limitations  are  inherent  in  the 
basic  assumptions  used  to  formulate  the  method.     Since  independent 
parameters  are  required,    a  great  number  of  practical  problems  can 
not  be  handled  by  this  method  as  it  is  presented  here.      By  extending 
the  mathematical  derivations  to  include  the  correlation  of  parameters, 
the  method  could  be  used  for  systems  in  which  the  parameters  are 
not  independent. 
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Table  I 
Functional  Formulas 


(1)  U  =  X  +  Y         py(u)  =  |px(x)pY(u-x)dx 

*  J°Px(u-y)pY(y)dy 

(2)  U  =  X  -  Y         Pu(u)  =  J~px(x)py(x-u)dx 

=  J,px(u+y)py(y)dy 

.  as 

(3)  U  =  XY  Pu(u)  =  J"px(x)py(u/x)^ 

-  •» 

(4)  U  =  X/Y  Pu(u)  =  -^2  Jpx(x)py(x/u)lxldx 

px(^y)py(y)iyidy 


0  pY(/u)  +   pv(-/u) 

(5)  U  =  X2  p..(u)  =  -* $ u^O 

U  2/u 


(6)  U  =  IXl  Pu(u)  =  px(x)  ♦  Px(-x)  x^O 


(7)  U  =  a/X  Pu(u)  =  -^  Px(a/u) 


(8)  U  =  aX  +  b        Pu(u)  =  -i-  px(^) 


(a  and  b  are  constants) 
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Table  II 
Sample  Data 


_^y 

41 

53 

65 

77 

89 

110 

0.9266 

0,7731 

0.6761 

0.6096 

0,5616 

120 

0.9204 

0.7669 

0.6699 

0.6034 

0.5554 

130 

0.9152 

0.7617 

0.6647 

0.5982 

0.5502 

140 

0.9108 

0.7573 

0.6603 

0.5938 

0.5458 

150 

0.9069 

0.7534 

0.6564 

0.5899 

0.5419 

4 


a 


34 
30 


41 


0.0165 


53 


0.00988 


65 


0.00658 


77 


0.00469 


89 


0.00352 


110 


0.000670 


120 


0.000561 


130 


0.000478 


140 


0.000413 


150 


0.000359 
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3 


Pu(u) 


U  =  XY 


Figure  1 


0.4- 


Pv(y) 


0.2- 


V  =  1/U 


10  11 


Figure  2 
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0.2" 


Pw(w) 


0.1" 


W  =   2V  +  1 


Figure  3 


0.2' 


PH(h) 


0.1- 


Figure  4 
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c 

START 

) 

N 

/ 

\       Read       / 

\   xit  p(xi),  a   / 

\    i=l,N     / 

\t 

i=l 

k=N 

U  =  a/X 


p(uk)=p(xi)xxf/a   ^ 


uk=a/xi 


i=i+l 
k=k-l 


Write 
uk,  p(uk) 

k=l,N 

\ 

1 

(  STOP   ) 

Calculation  of  u,  and  p(uk)  where  U  =  a/X  and  a,  x, 
and  p(xi)  are  known. 


Figure  5 
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(    START    ) 


i=i+l     6 


Yes 


CN   ? 
No 


(   STOP      ) 


Read 
i=l,N 


AX 


-,   Compute 
p(*J  by 
Interpolation 


p(x^)=p(xi) 


Write 
xr  p(Xj) 


xj+1  =  Xj  +  AX 


>  j-j+1 


To  perform  rectangular  integration  of  the  distribution  p(x. ) 
vs.  x. ,  the  distribution  must  be  approximated  by  rectangular 

This  program  calculates  the  x's 


areas  of  equal  widths  ax. 
such  that  x 


J+l 


=  x .  +  ax  and  interpolates  for  the  correspond- 
j 


ing  values  p(x.)  from  the  given  distribution  p(x.)  vs.  x., 


Figure  6 
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(    START   ) 


U  =  X  +  Y 


Read 
x±,   p(x±) ,   y^,    P(yj) 

Ax     (i=l.N:,1=l.M) 


l=i 


u(l)=x(l)-fy(I 


i=i+l 


u(i)=u(i-l)+Ax 


No 


Yes 


i=l 


p(u^)=0.0 


Yes 
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a  p(Uj)=p(Uj)+p(xi)xp(y^i+1)xA: 


Yes 


i=i+l 


Yes 


*  i=l 


Write 
uk,  p(uk) 

k=l,(N+M-l) 


Yes 


(  STOP  ) 


This  program  performs  the  convolution  U  =  X  +  Y  by  using 
graphical  convolution  techniques.   The  distributions  p(x.) 
vs.  x.  and  p(y^)  vs.  y.  are  represented  by  sets  of  rec- 
tangular  areas  of  width  ax.   The  distribution  of  x.  is 
reversed  and  passed  across  the  y  .  distribution  in  ax  in- 
crements.  The  distribution  p(u,)  vs.  u,  is  then  obtained 
by  the  normal  graphical  convolution  mathematical  calcula- 
tions [Bracewell  1965,  Healy  1969J  . 


Figure  7 
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C  START  ) 


U  =  XY 


i  =  l,n 
j  =  l,m 
k  =  l,b 


■*  k=k+l 


Read 
PY(yk)>  yk>  ^j 


Compute 

MS1 


I 


j=l 


6- 


Compute 


-1-1 


M-jH'd 


z=l 
k=l 


J 


pY(y;)=pT(yk) 


pY(yi)=o 


Compute 

pY(yi)    *>y 

Interpolation 


A  A 
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Write 

*r.<7;) 


Form 

PY(y^)] 


Yes 


No 


>  ^ 


z=z+l 


A   A 


Compute 


Compute 


Compute 

Pn(Uj)=TAX 


(  STOP  ) 


Figure  8 
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(  START  ) 


U  =  X/Y 


i  =  l,n 

j  =  l,m 
k  =  l,b 


*  k=k+l 


Read 
Py(yk),  yk,  u^ 


J=i 


5- 


i=l 


6* 


Compute 
*z   "  xi/uj 


Yes 


z=l 
k=l 


* 


pY(yi)=o 


^  i=i+l 


*pY(yi)=PY(yk) 


*  e* 
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± 

Compute 

Py(yi)    fey 

Interpolation 


Write 


e 


Form 


Yes 


^  z=z+l 


Compute 
-2 


u. 


Compute 
S  =  axu~2 


Compute 


Compute 
p^Uj^TS 


Compute 
T  =LR,PY(y^)] 


Write 
Pu(u.) 


»  d«j+i 


Yes 


jj 

(  STOP  ) 


Figure  9 
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APPENDIX  I 

The  solution  to  the  deterministic  problem  has  been  obtained  by- 
using  three  separate  programs.      Linear  interpolation  and  rectangular 
approximations  have  been  used  to  perform  all  convolution  and  rec- 
tangular integration. 

The  problem  and  given  information  is  as  follows: 
Obtain  the  cumulative  distribution  of  6  where 

9       =      C     +     X{*T    +    yfLj    4     wCw) 

C     =     0.175 

A     =     8.  1 

B     =     27.  8 

X(x)     is     N(130;10) 

Y(y)     is     N(65;12) 

W(w)  is  uniform  0.  875  from  -   4/7  to  +  4/7 

Program  one  inverts  the  distribution  x  vs.    p(x)  and  multiplies 
it  by  the  constant  A.      Program  two  does  the  same  for  B/Y(y).      Since 
both  X(x)  and  Y(y)  are  normal  distributions,    the  input  data  for  these 
programs  was  taken  from  standard  normal  tables.      25  data  points 
were  used  in  each,    the  median  and  each  1/4  standard  deviation  out  to 
3  standard  deviations.      The  resulting  distributions,    X  (x)  and  Y  (y), 

are  then  used  as  input  data  for  program  three.      Program  three  per- 

<  i  i 

forms  the  convolution  U(u)  =  X  (x)  +   Y  (y),    the  addition  U  (u)  =  U(u)  +  C, 
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and  the  convolution    9    =  U  (u)  +  W(w)  to  obtain  the  cumulative  dis- 
tribution of  0. 

To  obtain  highly  accurate   results  a  £x     of  0.  002  was  used  for 
convolution  and  rectangular  integration.     The  over-all  accuracy  for 
this  type  of  solution  will  depend  upon  how  closely  the  rectangular 
approximations   resemble  the  continuous  distributions  which  they 
represent. 
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ftPPPNOlX     T     -    PROGRAM    ONE 


C  rVMU^TION     IF     THE     FPCRA0ILITY    DENSITY     FUNCTION,  P(X«), 

C  WHERF     X»=A/X     \KC    A,     X,     AMP    P(X)     ARE    KNOWN.     ?  ( X I  HAS 

C  N(P>AL    DISTRIBUTION.     N=NUM8ER    OF    POINTS    X(I),     f*  IS    THE 

C  iTANDAPP    ['f-VIATION    OF    THE    DISTRIBUTION    P(X). 


100 
101 
102 

210 

211 

212 

213 

103 

214 

300 


301 
->  n  ^ 

201 


PI 
PL 
P  c 

Pt 

Ft 
FC 
FO 
WF 
FT 
WF 
FO 
WF 
FT 
WF 
FC 

on 

o  r 

Ff 

WP 

FO 

CT 

QO 

K* 

PL' 

U( 

cr 

WF 
FO 

PC 

FT 

ST 

r  r 


MP 

AF) 

AP 

AD 
RM 
P  M 

PM 

IT 
RM 
IT 
RM 
IT 
P  M 
IT 
PM 

"3 

AD 

FM 
IT 

FM 

NT 

3 
N- 
(K 
K) 
NT 
IT 
FM 

3 
RM 
OF 
P 


MSI 
(5, 

<5, 

(5, 

AT( 

AT( 

AT( 

E  it 

fill 

L  (t, 

AT( 

t  lb 

fill 

E  (6 

ATI 

00 

<c» 

AT( 

F  (6 

AT< 

1  Ml 

01 

1*1 

)  =  < 

=  A/ 
I  Mil 
F  (6 

AT( 

02 

AT( 


0*!    X(50) , PX(5C) ,U(5C) ,PU(50) 

1 C  0 )  N 

1  C  1  )  M. 

102)A 

12) 

IZ) 

F10.0) 

,210)     N 

1H1 , T20, 'N    =    '  ,  12//) 

,211  )    M 

T2C , 'M    =    • ,12//) 

,212)     \ 

T2C, 'A    =     «,F10.5//) 

,213) 

T20,«X(I)«,T5C,'PX(I)«//) 

1=1  ,N 

1  C  3  )     X  (  I  )  ,  P  X  (  I  ) 

2F1C .0) 

,21'  )     X( I )  ,PX( I  ) 

1 1 3  C  .  8  J 

F 

I  =  1  ,  N 

X(I )**2/A)*(PX( I )/M) 
X(I  ) 

f,  2C 0 1 

//,T2  0,'X«  ,T50, «P(X     )',//) 

1=1, N 

2F3^.B) 


N    =    2f 
M    =    10 

A     =     F.l 


PX 


C. 5  06  24 
^. 51428 
0.52258 
r  .  «=  3  1  1-+ 
C. 5  3  999 
C. 54Q15 
0.55862 
C. 54842 
7*57 
C. 58909 
C.c9999 
C. 61 132 
C.623C7 
C. 63529 
C. 64799 
0.66122 
C. 67499 
C.63c36 
0.70434 
C.71Q90 
0.73636 
0.75348 
0.77142 
0.79024 
C. 309^9 


C 


">Q63-01 
5t>3E-01 
05CE-01 
749E-01 
9^4~-01 
249E-01 
C62E-01 
100F-01 
137F-01 
085E-01 
995F-01 
070E-01 
636F-C1 
372E-01 
964E-01 
413E-01 
995E-01 
11CE-01 
749^-01 
967E-01 
353E-01 
7Q4F-G1 
835E-01 
374E-01 
970E-01 


0.  13 
0.27 
0.51 
0.91 
0.14 
C.23 
0.33 
0.45 
0.53 
0.70 
0.79 
0.83 
0.83 
0.77 
0.6  7 
0.55 
0.4  3 

0.31 
0.21 
C.  13 
O.HO 
0.45 
0.23 
0.11 
0.54 


3  0  6 1 6 
86»73 
90586 
01518 
99999 
17977 
61402 
77^79 
5  5  301 
27987 
22247 
81481 
22726 
60853 
92061 
782^7 
0222C 
12370 
14366 
4P<*36 
66663 
22629 
31943 
80331 
32096 


4E 
8E 
IE 
6E 
6E 
9E 
<5F 
4F 
4F 
7E 
3E 
9r 
4F 
6F 
7E 
IE 
2F 
3E 
lc 
3F 
7E 
7E 
7F 
2F 
7F 


01 
01 
01 
01 

02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
01 
01 
01 
01 
00 
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APPENDIX  I  -  PROGRAM  TWO 


C  EVALUATION  OF  THE  PROBABILITY  DENSITY  FUNCTION,  PIY'I, 

C  WNEPE  Y«=B/Y  AND  3,  Y,  AND  P(Y)  ARE  KNOWN.  ?(Y)  HAS 

C  NTPMAL  DISTRIBUTION,  N=NUMBER  OF  POINTS  Y(I),  M  IS  THE 

C  STANDARD  DEVIATION  OF  THE  DISTRIBUTION  P(Y). 

DIMENSION  X<50) ,PX(50) ,U<50) ,PU(50) 

RFAD(5,10C)N 

PtAD<5,101)M 

RtADI5tlC21A 

FOP  MA T( 12) 

Ff  -RM/T(I2I 

FORMAT  (Fl  CO) 

WRITE t 6, 2 JO)  N 

FCFMAT( 1H1 ,T20,'N  =  ',12//) 

WRITE(c,211)  M 

FORM/T(T20,,M  =  • ,12//) 

WRITE (6,212)  A 

F0RMAT(T2C, 'A  =  ',F10.5//) 


100 
101 
132 

210 

211 

212 

213 

103 

214 
300 


301 
200 
?01 


WRITER, 21 


) 


•Y(  I ) '  ,T50,  «PY(I  ) '//) 
N 
X( I ) ,PX(I ) 


FORMAT (T20, 

DO  300  1=1, 

PLAD(  5,103) 

RCFMAT(2F1C0) 

WPITE(5,21Z  )  X(I )  ,PX(  I) 

FORMAT<2E3C81 

CONTINUE 

DO  301  I=1,N 

K=N- I + 1 

PU(K)=(X(I)**2/A)*(PX(I)/M) 

U(K)=A/X< I ) 

CONTINUE 

WRITE (6,200) 

FCR^AT< //,T20, 

DO  302  1=1 tN 

FPFMAT(2E3C8) 

STOP 

END 


•Y»  ,T50, 'P(Y    )•,//) 


25 
12 
27.8 


Y« 

C.275 
0.283 
0.292 
C.302 
0.312 
0.323 
0.3  34 
0.347 
0.361 
0.375 
0.391 
0.^08 
0.427 
0.448 
0.471 
C.496 
C.524 
0.555 
0.591 
C.631 
C.678 
0.731 
0.794 
0.868 
0.958 


pyi 


24740 
67329 
63145 
173  73 
35939 
25566 
93960 
49985 
03874 
67550 
549U 
82331 
69212 
38685 
18622 
42831 
52803 
99976 
48908 
81788 
04  343 
57859 
28536 
74962 
62025 


00 
00 
00 
00 
00 
00 
0  0 
00 
CO 
CO 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 


0.13 

0.26 
0.47 
0.80 
0.12 
0.19 
0.26 
0.35 
0.43 
G.4° 
0.53 
0.53 
0.50 
0.44 
0.35 
0.23 
0.20 
0.  13 
0.85 
0.50 
0.27 
0.13 
0.64 
0.27 
0.11 


454551F 

19796<+E 

343379E 

423284E 

321751E 

1329  19E 

742373? 

031176E 

010139E 

425173E 

20552PE 

600159E 

520163E 

558611F 

740417E 

304B4*E 

377035E 

634053E 

751051E 

082946E 

210426E 

721460E 

261079F-01 

932845E-01 

092324F-01 


00 
00 
00 
00 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
00 
00 
00 
00 
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APPENDIX  I  -  PROGRAM  THREE 


CALCULATION  OF  U  =  X  +  Y  WHERE  X,  P(X).  Y,  P(Y),  ARE  GIVEN 
AND  THE  GIVEN  VALUES  CF  X(I)  AND  YU)  ARE  INCREMENTED 
PY  DX. 


131 
210 
1C6 
213 
105 
211 
104 
212 
202 

203 

302 


MENS  I  ON  X( 500) ,  PX ( 5CC )  ,  Y( 500 )  , PY ( 500  ) ,U ( 500 ) , PU ( 500  ) 
'(600) »Z(  1000) ,PZ(1000) ,c(500) ,W<600) 
A0(C,1C1)  OX 


DI. 

PK60C)  i  l\  rju'ji  t 

RF AD (5,101)     OX 

FCRMAT(FIO.O) 

WRITE(6t21C)     DX 

FCRM*T(1H1 ,T20, 'DX  =',F10.3//) 

RFAD( 5,106)  C 


FCRMAT(FIC.O) 
WF  ITE  U,213)    C 
rfPMAT(T20, 'C    =',F10. 
RtAD(5tl05)     N 


READ (5. 1051 

FCiRM^T(  15) 

Ml) 


B//) 


WFITF ( e ,21 1 )     N 

FCRM£T(T2C, »N    =•  , 15//) 

REA0(5,10«O     m 

FORMAT! 151 

WFITF (6,21 2)     M 

FOPMAT( //,T20, »V    =• ,15//) 

WRITE (6,202) 

FCPMAT(//,T20,'X(I),,T50,«PX(I)«//) 

CALL    UNF(X,PX) 
WPITE(6f203) 

FCPM/T(//,T2  0,,Y(I)',T50,,PY(I)«//) 
CALL    ONE(Y,PY) 


CALL     ONF(Y,PY) 
K  =  N=+M-1 

U( 1)=X( 1 )+Y( 1) 
nc    30  2    1=2, K 
U( I )=U( 1-1 )+DX 
DC    3C3    J=1,K 
PU( J  I =0.0 
DC    30^    1*1 t J 
I 


304 

303 


DC    30^    1=1,) 

IF( ( J-I+l) .GT.V)    GC    TG    304 

IF(  I.GT.N)    GO    TC    303 

PU(J)=PU(JH-PX(I)*PY<J-I  +  1) 

CONTINUE 


*ox 


CONTINUE 

CC  NTIMJF 

CALCULATION    of    L(J)=U(J)  +  C    WHEdE     »C»     IS    A    CONSTANT 


or   3C6   I = i . k 

IK  I  )=U( I )+C 
30  6    CONTINUE 

WRIT*-  (6.20CI 

200  FCRMAT(1H1,T20,«U(I),,T50,'PU(I)'//) 
nr  305  1=1, K 

305  WRITE(6.?01I  U( I ) ,PU( I ) 

201  FCP^T(2E3C8) 

CALCULATICK  AND  PLOT  CF  CUMULATIVE  DISTRIBUTION 

F(1)=PU( 1)*DX 

DO  320  1=2. K 

F(  I  )=F(  1-1  )  +  PII(  I  )  *0X 
3  20  CCNTINUF 

WPITF( ^,225) 
225  FCRMAT( 1H1 ,T20, 'U( J)  •  ,T40, 'CUMULATIVE  DISTRIBUTION*) 

DC  321  1=1 ,K 

WFITE(6,201 )  U( I ) ,F( I ) 
321  CCNTINUF 

PERFORM  THE  CONVOLUTION  OF  «U»  WITH  THE  CONSTANT  PDF 
•W«  (Z=U+W). 

CALL  TWC(W,PW,L) 
NM=K+L-1 
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310 


312 
311 

?30 

313 
231 


Z(  1)=U< 1)  +  K(1) 

or    310  1=2, NM 
Z<  I  )=Z( 1-1 )+DX 
or  31  1  J=l ,NM 
P7< J)=O.G 
DC  312  1=1  ,  J 

if( ( j-T+i) .gt.d  go  rn  312 

IF(  I.CT.K  )  GO  TC  311 

F7( J)=P7( J)+PU( I )*PW( J-I+1)*0X 

CONTINUE 

CONTINUE 

WRITE (6, 230) 

FCPMAK LH1  ,T20,  «Z( I ) •  ,T50,'PZ(I  )  »//) 

00  313  1=1 ,NM 

WPITFU>,2?1)  Z ( I )  ,  PZ (  I  ) 

FCRMATt2E30.81 

CALCULATION  OF  THE  CUMULATIVE  DISTRIBUTION  OF  »Z« 


PZ(  1)=PZ(1 )*DX 
00  33C  1=2, NM 
PZ( I )=FZ( I-1)+PZ( I )*DX 
330  CONTINUE 

WFITF(6,220) 

220  FDPMAT(1H1,T20, »U(I  )»  ,T*0, 'CUMULATIVE  DISTRIBUTION*) 
DO  331  T=l ,NM 

221  FCRN,AT(2E3C3) 
3  31  CONTINUF 

cTOP 
EMO 


sueeouTir-t   on£(xe,pxej 

FROM  GIVEN  VALUES  OF  X(I)  AND  PX(I),(THESE  ARE  U(I) 
AND  PU(I)  PREVIOUSLY  CALCULATED),  AND  A  CHOSEN  OX,  WE 
OBTAIN  THF  VALUES  OF  X(J),  WITH  INT5VALS  CX,  AND 
C;LCULATE  THE  CORRESPONDING  VALUES  OF  °X(J)  by  linear 
INTERPOLATION* 

DIMENSION  X{500) ,PX(5C0) ,XE(SOO)  ,PXE(500) 
RE AG< 5,1001  N 

IOC  PC'RMAT(  12) 

PTAD(r, 101)  OX 
101  FC PMAT(FIO.O) 

r\n   inn   T  _  l   m 

300 
102 


400 
502 


501 

5C0 

2C0 
201 

203 


DO  300  1=1. N 


UU  5UV     1  =1  ,  IN 

P  F  A 0  (  5  ,  1  C ?  )  X( 

P0RMAT(2F10.0) 
1  =  1 


<  I  ),  PX( I ) 


J=l 

XF  (  J)=X( T )+0X/2 

IP(XE( J)-X(I  )  )  500,5C1,5C2 

1  =  1  +  1 

IF( I.LT. (N+l) )  GO  TC  400 

GO  TO  203 

PXE< J)=PX( I ) 

GO  TC  200 

PXE( 


PXE( J)=PX( I  ) 

GO  TC  200 

PXE(J)=FX(I-1)+(PX(I)-PX(I-1))/(X(I>-X(I-1))*(XE(J>- 

X(  !-l)  ) 

WRITE(6,20IJ   XE(J),  PXE(J) 

FORMAT! 2E30. 81 

XF( J+1)=XE( J)+DX 


XF( J+1)=XE( J)+DX 
J  =  J+1 

GO  TC  ^00 
RETURN 

END 


SUBROUTINE  TWO(W,PW,L) 

C      THIS  PROGRAM  TAKES  THE  GIVEN 
C      FUNCTION  'W«,  WITH  AMPLITUDE 


CONSTANT  PROBABILITY 
0.P75  FROM  -4/7  TO  +4/7, 


43 


AND  CPTAINS  THE  VALUES  GF  W(I)  WITH  INCREMENT  CX=0.C02 


25C 


3  50 


251 
351 


DIMENSION  W(600) 

W(l)=-0.571 

FK1  >=0.875 

PX =0.00  2 

L  =  571 

WRITE (6, 250 ) 

FfPMAT( 1H1 , T20,  • 

PT  ?50  1=2, L 

W(  T  )=W(  I-D+9X 

FW(I  )=PW( 1-1) 

CONTINUE 

DP  3el  1=1 ,L 

WF ITE (6,251 )  W( I ) 

FC-RMAT(2E30.  <?) 

CC  NTTNUF 

RETURN 

f  MP 


PW(600) 


W(  I  )  •  ,T50»«PW(  IP//) 


PMI  ) 


N  =  15 

C  =  6.175 

DX  =  C.002 
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THETA 


CUMULATIVE 
CISTRIPUTIOM 


•0,63128 

•C.65128 
■C. 64128 
■C.  62128 
•C. 60128 
■C. 58128 
•0.56128 
■0.54128 
■C. 52128 
■C. 50128 
■C. 48128 
■C. 46128 
■0.44128 
■0.42128 
■C.40128 
■0.33123 
36128 
34  128 
32123 
30128 
23128 
26128 
2^128 
22123 
20129 
18128 


-C. 16128 


14126 
0.12128 
0.10128 
C.  812*30 

C. 61280 
C. 41290 
0.21280 
0.12303 
0.13719 
C. 38719 
0.53719 
C.73719 
0. 9  8719 
C. 11871 
0.  13871 
C. 15  871 
C. 17871 
0.19871 
0.21671 
0.23871 
C.25S71 
C27871 
0.29871 
C. 31371 
0.33971 
0.35871 
0.37871 
0.39871 
0.41871 
0.43871 
0.45871 
0.47871 
0.49971 
0.51871 
0.53871 
0.55871 
C. 57871 
0.59671 
0.61871 
0.63871 
0.65871 
0.67871 


109E-01 
075E-01 
0416-01 
045E-01 
044F-01 
0446-01 
044E-01 
343P-01 
0436-01 
043--01 
042E-01 
042E-01 
042E-01 
041F-Q1 
0416-01 
041E-01 
0406-01 
040F-01 
040E-01 
0^96-01 
039E-01 
039E-01 
0386-01 
0336-01 
0386-01 
0376-01 
03  76-01 
0376-01 
0366-01 
0366-01 
3536-02 
355E-C2 
3526-02 
3536-02 

644E-02 
6436-02 
628E-02 
5945-02 
5606-02 
3536-01 
3496-01 
Q46F-01 
9426-01 
939F-01 
9366-01 
9326-01 
^ 96-01 
9256-01 
9226-01 
9196-01 
9156-01 
9125-01 
908F-01 
9056-01 
9C2F-01 
9986-01 
8956-01 
891E-01 
968E-01 
885E-01 
3815-01 
978E-01 
8746-01 
3716-01 
368E-01 
3615-01 
3355-01 
8096-01 


0.351622805-08 
0.2^3545596-07 
0.933422416-07 
0.293669135-06 
0.71189C306-06 
0.  1*7910676-05 
0.273452416-05 
0.46203268E-05 
0.72799021E-05 
0.10959O82E-04 
0. 155095296-04 
0.213390276-04 
0.236690055-04 
0. 375373066-04 
0.4319f212E-04 
0.603560955-04 
0.75755*956-04 
0.931420656-04 
C. 113283376-03 
0.  13646450E-03 
0.  162986976-03 
19316789E-03 
227340216-03 
265353016-03 
309071276-03 
357373626-03 
0. 411151216-03 
0.47080661F-03 
0.536753336-03 
0.60941558E-03 
0.699226675-03 
0.776627306-03 
0.872066015-03 
0.975992535-03 
0.103886166-02 
0.  121112946-02 
0.134325356-02 
0.  148569236-02 
0.16399030E-02 
0.130333906-02 
0.  197944776-02 
0.216766876-02 
0.236843416-02 
0.258216846-02 
0.230928916-02 
0.30502048E-02 
0.330531226-02 
0.357499286-02 
0.38594098E-02 
0.415950646-02 
0.447501995-02 
0.430646646-02 
0.51541552E-02 
0.551338066-02 
0.539941446-02 
0.62975064E-02 
0.67128837E-02 
0. 714575136-02 
0.75962991E-02 
0.306470966-02 
0.855114316-02 
0.905575235-02 
0.957866395-02 
0.101199866-01 
0.106797965-01 
0.112581516-01 
0.118550926-01 
0.124706515-01 
0.  13104349E-01 
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C.69P71783E-01 

0.71871753E-01 

0.73971732E-01 

C.758717C6E-C1 

C.77871680E-01 

C.79371655E-01 

C.31871629E-01 

C. 33871603E-01 

C.85&71577E-01 

C.«7671552F-01 

C.89971526E-01 

r.<U8715C0E-Cl 

C.93871474E-0I 

C.95671449E-01 

0.97871423F-01 

f.99«71397E-01 

C.10187137E  OC 

C.10387135E  00 

C.1D537132E  00 

C.10787129E  CO 

C.10987127E  00 

r,11187124E  00 

0.11387122E  00 

C.11987119E  CC 

C.11787117E  00 


0.  11997114E 


10 


C.12137111E  00 
C.123971C9E  00 
C.125371C6E  00 
G.127871C4E  00 
C.129371C1E  OC 
G.13187C99E  00 
0.13337096F  00 
C.13587093E  00 
G.13787091E  OC 
C.13987083E  CO 
C.14197336E  00 
C.14387083E  CO 
C.14587090E  00 
C.14797078E  OC 
C.14987075E  00 
C.15197C73E  00 
0.15387070E  CO 
G.ic;c87068E  GO 
C.15737C6EE  OC 
C.15Q87C62E  00 
C.1618706CE  00 
0.16397057E  00 
0.16537055E  CO 
C.16737052f  CO 
T.1698705CE  00 
C.171870^7F  OC 
C.17387044E  00 
r.l75«7042E  00 
C. 17787C39E  00 
0.17Q37037E  00 
C.13187034E  00 
C.13387032E  CC 
C.13587029E  00 
0.137R7026E  00 
0.18967024L  00 
0.19187021E  OC 
C.19367019E  00 
C.19567016E  00 
C.19787014E  CO 
C.19937011E  00 
C.20187008E  OC 
C.203870C6E  03 
C.205870C3E  00 
0.20787001E  00 
C.20986998E  00 
C.211»6996c  00 


0.13757691E-01 

0.14429174E-01 

0.  15119277E-01 

0.15Q279E0E-01 

0.  1S555130E-01 

0.17300732E-01 

0.  13064659E-01 

0.19846810E-01 

0,  19647073E-01 

0.20465329E-01 

0.21301460E-01 

0.22155326E-01 

0.23026783E-01 

0.23915663E-01 

0.2^9217926-01 

0.2574A982E-01 

0.26665044E-01 

0.27641781F-01 

0.28614990E-01 

G.29604476E-01 

0.33f>10032E-01 

0.31631447E-01 

0.32668509F-01 

0.337209O6F-01 

0.34798663^-01 

0. 35871336E-01 

0.3A968727E-01 

0. 33080610E-01 

0. 29206754F-01 

0.40346917E-01 

0.A1500859E-01 

0.4?S63343E-r»l 

0.43349129F-01 

0.45042980E-01 

0.AS2^C?662E-01 

0.47468934E-01 

0.A9700567E-01 

0.49944323E-01 

0.51199969F-01 

O.52467272E-01 

0.53745996E-01 

0.55035912E-01 

0.5S3367P3E-01 

0.57*48379E-01 

0.58970470E-01 

0.60302B28E-01 

0.61645232E-01 

0.62997460E-01 

0.64359307E-01 

0.65730572E-01 

0.67111015E-01 

0.63500519E-01 

0.69898844E-01 

0.71305752E-01 

0.72721124E-01 

0.74144721E-01 

0.75576365E-01 

0.77015817E-01 

0.734623C9E-01 

0.79917431E-01 

0.81379235E-01 

0.fi2348191E-01 

0.84324062E-01 

0.85306727E-01 

0.87296009E-01 

0.89791728E-01 

0.902°?765E-01 

0.91802^01E-01 

0.93316257E-01 

0.9A93^414E-01 

0.9S362293E-01 

0.Q7393715E-01 
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0.213*69935  00 

0.21586990E  00 

C.217P6988E  00 

0.21 986 985 E  00 

0.22186583E  00 

0.223S6980E  00 

C.22586977E  00 

0.22766975c  00 

C.22<=86972t  00 

0.23136970E  00 

G.23336967E  00 

C.23586965c  00 

0.237869625  00 

C.23(:i8695::>E  00 

0.24186957E  00 

0.24386954E  00 

0.24586952E  00 

0. 24786 949£  0  0 

C.24986947E  00 

C.25  186944E  00 

C.25386941E  00 

C.25C<J6^39E  00 

C.25"786936E  00 

0.25986934E  03 

0.261 86 931 E  00 

G.26386929E  30 

C. 265869265  00 

C.26786923E  00 

C.26986921c  CO 

C.27186913E  00 

C.27386916E  00 

C. 275869135  00 

C. 27786 91  IE  00 

C.279869C8E  00 

0.23186905E  00 

C.283869C35  OC 

0.28586900F  00 

0.28786898E  00 

0.28986895E  00 

0.29186893E  OC 

0.29386890E  CO 

0.29586887E  OC 

0.297868855  03 

0.29986882E  00 

0.30186880E  00 

C.30386877E  00 

C.30536374E  OC 

0.30786372E  00 

C. 309868695  00 

0,311668675  00 

C. 313868645  00 

C.31586862E  00 

0.31786859E  OC 

C.31936856F  00 

0.32166354E  00 

0.323868515  00 

0. 32586649E  00 

0.32786346F  00 

0.32986344E  00 

0.331868415  00 

0.33386838E  00 

0.33586836E  00 

0.33786833E  00 

0.339863315  00 

0.34186328c  00 

0.34366826E  00 

C.345P6823E  00 

0.3478682CE  00 

0.34986818E  00 

C.35186315E  OC 

0.353868135  OC 

0.3558631CE  00 


99<V30561E-01 


10097271E 
102519C9E 
13407227E 
10562944F 
1 3  7  1 9 1  3  8  F 
13375791E 
11032397E 
1119C450E 
11343432F 
11506838E 
11665654E 
11824870E 
11984479E 
12144470F 
12304831E 
12465549E 
126266185 
12738033E 
12^497775 
13111842E 
13274223E 
13436908F 
13599885F 
13763154E 
13926703E 
14090532E 
14254630E 
14418989E 
14883606E 
1474R472E 
14913583E 
15078°32E 
152445145 
15410322F 
15576351E 
157425945 
15309046E 
160757015 
16242546F 
16409582E 
16576803E 
16744202E 
16911775E 
17079520E 
17247432E 
17415506E 
17583734F 
17752123E 
17920661F 
18089348E 
132581785 
13*271515 
13596262F 
18765509F 
18334888E 
19104391E 
19274020E 
13443768E 
19613636F 
19783616E 
19953704E 
20123899E 
202942C1F 
234645995 
23635098E 
20305693c 
23976377E 
21147156E 
213180185 
21438971E 
21660006E 


00 
00 
00 
00 
OC 
00 
00 
CO 
00 
30 
00 
00 
00 
00 
00 
00 
00 
OC 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
OC 
00 
00 
00 
00 
00 
00 
00 
00 
OC 
00 
00 
00 
00 
00 
00 
00 
0  0 
00 
00 
00 
0  0 
00 
00 
00 
30 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 


47 


c 
c 
c 
c 
r 
r 
C 
0 

c 

r 

C 

c 

o 
o 

r, 
0 

c 

0 

c 
c 

Q 

o! 
c. 

c. 

r, 

c 

c 

c, 

0, 

c. 
c 

0. 

0, 
0, 

c\ 
c, 
c 

c, 


357 
35^ 
361 
363 
365 
367 
369 
371 
373 
375 
377 
379 
381 
383 
335 
33  7 
339 
391 
3^3 
395 

3  97 
39° 
-01 

,403 

,405 

,407 

,40G 

411 

413 

415 

,417 

,41  Q 

421 
42  3 
425 

427 

u.7C. 

431 

433 
435 
437 
439 
441 
44? 
4'.  5 
447 
h4<~- 
451 

453 
455 
457 

461 
46? 

4  65 
67 
6C 
71 
73 
75 
77 
70 

1 

3 

c. 
7 
c 


C.49 
C.49 
0.49 
C.49 


868C3 
86805 
368C2 
86300 
36797 
36795 
86792 
36790 
86787 
86784 
86^62 
86779 
86777 
66774 
06771 
86769 
86766 
86^64 
86761 
86759 
86756 
c6753 
P6751 
86743 
P67n6 
66  74  3 
86741 
86  7  33 
86735 
86733 
8673C 
86723 
8fe7?5 
86723 
P672C 
86717 
36715 
36712 
86710 
»67C7 
8b7C5 

867C2 
8  6  6  9  9 
86697 
86694 
0  6  6  9  2 
86689 
85687 
36684 
86681 
8667^ 
«6676 
86674 
86671 
°6668 
86666 
86663 
86661 
86658 
86656 
66653 
8665C 
86o43 
86545 
86643 
866^  C 
P6o38 
«6635 
«6632 
36S3C 
86627 
86625 


6  00 
c  00 
E  00 
E  00 
E  00 
F  00 
F  00 
E  00 
E  00 
t  00 
5  00 
F  00 
E  00 
E  CO 
E  00 
E  CO 
E  00 
c  00 
E  00 
F  00 
E  00 
F  00 
E  00 
E  OC 
F  00 

E  00 
00 
CO 
00 
00 
00 
CC 
CO 
00 
CO 
00 
00 
CO 
00 
OC 
CO 
0  0 

E  00 

E  00 

E 

E 

E 

E 

c 
f: 
E 


00 
00 
CO 
CO 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
OC 
00 
00 
00 
00 
CO 
00 
00 


0. 

21331125E 

OC 

0. 

22002327E 

00 

0. 

22173607E 

00 

c. 

22344971F 

00 

0. 

22516406F 

OC 

0. 

22687916E 

OC 

n, 

228595C2E 

00 

0. 

23031157E 

00 

0. 

23202884E 

OC 

0. 

23374677E 

00 

0. 

23546535E 

00 

0, 

23713458E 

00 

0. 

23390442P 

00 

0, 

24062A85E 

OC 

0. 

24234587E 

00 

P. 

2*4C6749F 

00 

0, 

24576Q65E 

OC 

0. 

24751234E 

00 

0, 

2*t923557F 

0  0 

0, 

25095928E 

00 

0, 

25268352F 

00 

0. 

25440824F 

OC 

0. 

25613338E 

00 

0, 

25785899E 

00 

c. 

2595«5C8E 

00 

c. 

251311591: 

OC 

0. 

2S303852E 

00 

n. 

26476586E 

00 

c . 

25549362F 

00 

0, 

.26622180E 

00 

0. 

26995039E 

00 

0. 

27167934E 

00 

0. 

27340871E 

00 

0. 

?75138^4E 

00 

0. 

27686852E 

00 

0. 

27359396E 

00 

r. 

23032976E 

OC 

0, 

23206092E 

00 

0  • 

28379244F 

00 

0, 

23552425E 

00 

0. 

23725636E 

00 

0. 

233c;88e3E 

00 

0. 

29072160E 

00 

0. 

29245460E 

00 

0. 

29413790E 

00 

0, 

29592150E 

00 

r . 

29765534E 

00 

~. 

29938948E 

OC 

0. 

30L12386E 

00 

0, 

30285847E 

0  0 

0. 

3045C332E 

00 

0. 

30632842E 

00 

0, 

303O't375E 

00 

0, 

30979931E 

00 

0, 

3U53506F 

00 

0, 

31327105F 

00 

c. 

31500721F 

00 

0, 

31674361P 

00 

0. 

3134801QE 

00 

0. 

32021695E 

00 

0. 

321953Q5F 

00 

0. 

32369113F 

00 

0. 

32542849F 

00 

0. 

32716602E 

00 

0. 

32390373F 

00 

0. 

33064163E 

00 

0. 

33237970F 

00 

0, 

33411795F 
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c. 

c. 
c. 

0. 
0. 

c. 

0. 

c. 
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c. 

0. 
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c. 
c. 

0. 

0. 
0. 
0. 

c. 

0. 
0. 

0. 

c. 

c. 

0. 
0. 

0. 

c. 


077 
079 
.081 
0  33 
085 
.0  87 
.089 
091 
0  93 
095 
097 
099 
101 
103 
105 
1C7 
109 
111 
113 
115 
117 
119 
121 
123 
125 
127 
129 
131 
133 
135 
137 
139 
141 

1*3 

145 

147 
149 
151 
153 
155 
157 
159 
161 
163 
165 
167 
169 
171 
173 
175 
177 
179 
181 
183 
185 
187 
189 
191 
193 
195 
197 
199 
201 
203 
2C5 
.2  07 
.209 
211 
213 
215 
217 
219 


8542E 
8540E 
8539E 
S537E 
6536E 
8534E 
8533E 
8531E 
8530E 
8529E 
8527F 
8526E 
8524E 
85  23E 
8521E 
8520E 
8518E 
3517E 
8515c 
8514F 
8513E 
8  5 1 1 F 
8510c 
3508E 
8507E 
8  5  05E 
8504E 
8502E 
8501E 
6500E 
8498E 
3497E 
8495E 
849^E 
3492F 
8491E 
8489E 
8488F 
8486E 
3485E 
8484E 
84826 
8481E 
8479E 
8478E 
3476E 
8475E 
8473F 
8472E 
8471E 
8469F 
8468E 
8466E 
8465F 
8463F 
8462  £ 
8460F 
8459E 
8457E 
8456E 
8455E 
8453E 
84525 
8450F 
8449E 
8447E 
3446E 
84<^E 
8443E 
8442E 
8440E 
3439F 


CI 
01 

01 
01 
CI 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
CI 
01 
01 
01 
01 
01 
CI 
01 
01 
01 
01 
01 
01 
CI 
CI 
01 
01 
01 
01 
01 
CI 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
CI 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
CI 


C.84500587E  00 
0.84674984E  00 
0.84849358E  00 
0.S5023695F  00 
0.85197985F  00 
0.85372216E  00 
0.85546368F  00 
0.85720426E  00 
0.85R94376F  00 
0.86C68207F  00 
C.8S241895E  00 
0.86415^22E  00 
0.8653877GE  00 
0.86761922F  0C 
0.36934847F  00 
0.87107527E  OC 
0.37279928F  00 
0.87452024F  00 
0.67623787F  00 
0.37795186F  00 
0.87966186E  00 
0.83136750E  00 
0.33306844E  00 
0.88476431E  OC 
0.33645470E  00 
0.88813919F  00 
0. 8  3  931 742 E  00 
0.89148391E  OC 
0.39315325F  OC 
0.39481002F  00 
0.39645375F  00 
0.89809901E  00 
C.39973032E  00 
0.93135223E  00 
0.90296429F  00 
0.91456605F  00 
0.906157C2E  00 
0.90773672E  00 
0.90930480E  00 
0.910S6072F  00 
0.91240412E  00 
0.91393453F  00 
0.91545159F  00 
0.91695482F  00 
0.91844386E  00 
0.91991836E  00 
C.92137790E  00 
0.92282218E  00 
0.92425084E  00 
0.92566353E  OC 
0.92706001E  00 
0.92843997F  00 
0.92980313E  00 
0.9311^918E  00 
0.932477Q5E  00 
0.93378925F  00 
0.93508285E  00 
0.C53535363F  00 
0,937fcl635E  OC 
0.93385589E  00 
0.94007713E  00 
0.9U27995E  00 
0.94246429F  00 
0.94363010E  00 
0.9<t477731E  00 
0.94590592F  00 
0.94701588E  00 
0.94310718E  00 
0.94917989E  00 
0.95023394F  00 
0.95126O45E  00 
0.95223642E  00 
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r, 
0. 

c, 
o, 

0, 
0, 
0, 

0, 

r 

0, 
5, 

0, 

c 

0, 
0, 


221 
223 
.225 
227 
229 
231 
233 
235 
237 
239 
2M 
243 
245 
247 
2^9 
251 
253 
.255 
.257 
259 
261 
263 
.265 
267 
269 
271 
273 
275 
277 
279 
281 
263 
285 
267 
269 
291 
293 
295 
2^7 
299 
3C1 
303 
3C5 
307 
309 
311 
313 
315 
317 
319 
321 
323 
325 
327 
3  29 
331 
333 
335 
337 
339 
341 
3*3 
345 
347 
349 
351 
353 
355 
357 
359 
361 
363 


9437E 
6436E 
e434E 
8433c 
8431F 
8430F 
842*5 
34276 
8426E 
8424E 
8423E 
8421F 
842CF 
8413E 
8417F 
841  5  F 
3414E 
8413F 
8411E 
8410E 
34085 
8407F 
3405E 
3404E 
3402E 
8<fClE 
34C0E 
8398F 
8397F 
3395E 
8394F 
8392F 
3391E 
338QE 
«383F 
8386E 
8385E 
6334c 
6  3  82E 
3381E 
6379E 
8373E 
3376F 
3375F 
8373E 
8  3  72F 
8371F 
3369F 
8363'= 
836bF 
8365E 
8^63E 
3362E 
8360E 
8359F 
»357E 
R356F 
8355E 
3353E 
8352E 
8350E 
8349F 
8347E 
3346E 
3344E 
3343E 
8342F 
S340E 
8339E 
8337F 
8336F 
8334E 


CI 
CI 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
CI 
CI 
01 
01 
01 
CI 
CI 
01 
01 
01 
01 
01 
01 
CI 
CI 
CI 
CI 
01 
01 
CI 
CI 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
CI 
01 
01 
CI 
CI 
CI 
CI 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 


0 

.95328496F 

00 

0, 

.95426524F 

00 

0, 

.95522726F 

00 

0, 

.95617115E 

00 

0. 

.95709705E 

00 

0, 

.95800507E 

00 

0. 

.95389539E 

00 

0, 

.95976812E 

00 

0, 

•96062338F 

00 

0, 

,961461*2F 

00 

c, 

►96228242F 

00 

c. 

.96308655E 

00 

0, 

,963873985 

00 

0, 

.954644C1F 

00 

0. 

,96539956c 

00 

0, 

•96613818P 

00 

0, 

.96636095E 

00 

0. 

.96756804F 

00 

0, 

.96325969F 

oc 

0, 

►96393615E 

00 

0. 

,9695976^F 

oc 

0. 

,970244416 

00 

0, 

.97097669E 

00 

0, 

.97149473E 

oc 

0. 

.97209877F 

00 

0, 

•97263397E 

00 

c, 

,97326«565E 

00 

0, 

,97382^03F 

00 

0. 

,9743^0366 

00 

0. 

.97491638F 

oc 

0. 

,9754<*175E 

00 

0, 

.97595429E 

00 

0, 

,976t5468F 

oc 

0, 

•97694319E 

00 

c« 

.97742003E 

00 

r. 

,977995485 

00 

o, 

.97833973F 

00 

0, 

.97378301E 

00 

0, 

.97921556E 

00 

0, 

.97963762F 

00 

0, 

,930049435 

00 

o, 

.93045117P 

00 

0, 

,930943075 

00 

0, 

.93122537F 

00 

0, 

,Q8159826E 

00 

0, 

,9819ei91F 

00 

c. 

.98231655F 

oc 

0. 

.93266238E 

00 

0, 

.93299-156E 

00 

0, 

,°9332834E 

00 

0, 

,9d364890F 

00 

0. 

.93396140E 

00 

0, 

.9342661CF 

oc 

c. 

.93456317E 

00 

0. 

93485273F 

00 

J. 

.93513502E 

00 

c, 

93541021F 

00 

0« 

,985678495 

00 

0, 

93593998F 

00 

0, 

93619485F 

00 

•98644322F 

00 

0. 

93663528F 

00 

0. 

98692119E 

00 

0, 

93715109E 

00 

0, 

9873750eF 

00 

0, 

93759335E 

00 

0. 

.937S0602E 

00 

0, 

93:<01327E 

00 

c. 

93821527F 

00 

0. 

93841208F 

00 

0. 

93960395E 

00 

0. 

93379093F 

00 
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365 

367 

369 

371 
3  73 
375 

377 
37o 

381 
383 
385 
387 
389 
391 
393 
395 

397 
30c 

401 

403 

^05 

407 

.409 

.411 

.413 

415 
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U  3.0 

.421 

.423 

.425 

.427 
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433 
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,437 

.439 

.441 

.443 

.445 

.447 
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.461 

.463 
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.469 

471 
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475 

.477 

479 

.481 

483 
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487 
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.4^1 

4C3 

!497 
.499 
501 
503 
505 
507 


3333E 
3331E 
333CE 
8328E 
8327E 
8326E 
8324E 
8323E 
8^21E 
332CE 
8318E 
8317E 
8315E 
8314E 
8313E 
8311C 
331CC 
8308F 
3307E 
83  05E 
83C4F 
3302E 
8301F 
8299E 
329PC 
3297E 
3295E 
3294E 
8292F 
8291E 
3289E 
3288E 
8286E 
8285F 
3284E 
8282E 
8281E 
8279E 
8278E 
8276c 
8275E 
8273E 
8272E 
8270E 
8269E 
8268E 
3266F 
8265E 
8263E 
3262F 
8260E 
8259E 
8257E 
8256E 
3255F 
3253E 
3  2  52E 
8250F 
8249E 
8247E 
8246E 
8244E 
8243F 
3242c 
8240E 
3239E 
3237E 
8236F 
3234E 
8233E 
8231E 
3230E 
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01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
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01 
01 
01 
01 
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01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
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01 
01 
01 
01 
01 
01 
01 
01 
01 
CI 
01 
01 
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01 
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01 
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0.9 
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0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
G.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
C.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
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0.9 
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0.9 
0.9 
0.9 


3897314E 
3915070E 
3932374E 
3949236E 
3965669F 
8981A76F 
3997277E 
9012476E 
9027282E 
9041706E 
9955761E 
90&9452E 
9082792F 
9095792E 
910Q<r64E 
9120314F 
9132846E 
914457SF 
9156016E 
9167162F 
9179028F 
9188614E 
9198C31F 
9208987E 
9213786E 
9223334E 
9237639F 
9246705E 
9255538E 
9264145E 
9272531E 
9280703E 
9238666E 
9296427E 
9303997F 
9311376F 
9318570E 
9325585E 
9332426E 
9339096E 
9345601E 
9351943E 
9358124E 
9364150E 
9370027E 
9375755E 
9381340E 
9386781E 
9392086E 
9397254E 
9402291E 
9407202E 
94119P8F 
9416655E 
9421203E 
=>425633F 
9429965E 
9434185E 
9438298E 
9442309E 
9446225F 
9450046E 
9453777F 
9457413E 
9460959E 
9464422E 
9467796E 
9471086F 
9474299E 
9477428E 
9480430E 
9483454E 
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C.15C98228F  01 
C.15118227E  01 
C.15138226E  01 
0.151^8224F  01 
6.15173223F  01 
C.151C3221F  01 
C.15218220F  01 
0.1523821PF  01 
0.152^8217E  01 
0.15273215E  01 
C.15298214E  01 
C.15318213E  01 
C.15338211F  01 
C.15353210E  01 
0.15378208F  31 
C.15398207F  01 
C.154182C5E  01 
0.15438204-E  01 
0.15^532029  01 
C.15478201E  01 
0.1549Q199F  01 
0.15518198^  01 
C. 1C538197E  01 
C.15558195F  01 
0.15578194?.  01 
0.155981 92 E  01 
0.1^618191F  01 
C.1563318CF  01 
C.15658188F  01 
C.15678186E  01 
C.15t9818^E  CI 
0.15718184E  01 
0.15738182E  CI 
0.15758181?  01 
0.15778179E  01 
0.15798173E  01 
C.15618176c  01 
0.1583S175E  01 
C.15858173F  CI 
0.15378172E  01 
C.15898170E  01 
C.1591B169F  01 
C.15933163E  01 
0.15958166E  01 
0.15978165F  01 
0.15C98163E  01 
C.16018162E  01 
G.16C38160E  01 
C.15C53159F  01 
0.16076157E  01 
C.16098156F  01 
0  .  1 5 1 1 8 1 5  5  E  01 
C.16138153E  CI 
C . 1 5 1 5  6 1 5  2  E  01 
C.15173150F  01 
0.16198149E  01 
C.16218147E  01 
0.16238146E  01 
C1625S144?  01 
C.16278143F  01 
C.16298141E  01 
P.16313140F  01 
T.16338139E  CI 
C. 16358137c  01 
0.16378136E  CI 
0.16393134E  CI 
C.16413133F  01 
C.16438131F  01 
0.16458130'  01 
0.164791286  01 
C.164C8127F  01 
C.15513126E  01 


0.99486357E  00 
0.99499188F  00 
0.99491Q46F  00 
0.99494636F  00 
0.99497259E  00 
C. 994998169  00 
0.995023C7E  00 
0.93504739E  00 
C.99507105E  00 
0.99509A18F  00 
0.99511671F  00 
0.99513870F  00 
0.99516016F  00 
0.99518108E  0C 
0.99520153E  00 
0.99522144F  00 
C.99524087F  00 
0.99525982?  00 
0.99527836F  00 
0.99529642E  00 
0.99531406c  00 
0. 9953312 9E  0  0 
0.9953^<U0E  00 
0.99536449E  00 
0.99538046E  00 
0.995396C2F  00 
0.93541116F  00 
0.99542594E  00 
0.9954^0369  00 
0.93545^43F  00 
C.995*6;>14F  00 
0.995481^99  00 
0.995H9448E  00 
0.99550712E  00 
0.99551946F  00 
0.99553150F  00 
0.99554318E  00 
D.99555457E  00 
0.99555565F  00 
C.9955765CE  00 
0.995537C5E  00 
0.99559730E  00 
0.9956C732F  00 
0.99561709E  00 
0.99562663F  00 
0.99563593E  00 
0.99564499E  00 
0.99565381E  00 
0.93566239E  00 
0.99567080E  00 
0.9:'5578o6E  00 
0.99568695F  00 
0.995^9470E  00 
0.99570227F  00 
0.99570960F  00 
0.99571675E  00 
0.99572372F  00 
0.99573052E  00 
0.C9573714E  00 
0.99574357F  00 
0.99574':>e3F  00 
0.99575591F  00 
0.99576181F  00 
0.99576753F  00 
0.99^77308?  0^ 
0.9?577844E  00 
0.99578369c  00 
0.99578*759  00 
0.99579370F  00 
0.99579347E  00 
0.99580312c  00 
0.99580759F  00 
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c. 
c, 

c\ 
c. 
c 


c 

c. 
c, 


C.16533124F  01 
0.16^81235  01 
C.16578121E  01 
0.1659812GE  01 
C. 166131185  CI 
C.16633117F  CI 
0;16658115E  01 
C.16678114E  CI 
0.16693112E  01 
0.16718111E  01 
0,167381106  CI 
0.167581C8E  CI 
0.1677P1C75  01 
C.167981C5E  01 
0.168131C^t  01 
C.16P331C2E  CI 
C.16858101E  CI 
L6878Q99E  01 
16898098E  CI 
L6918097F  01 
169330956  01 
16958094E  01 
C.16978092E  CI 
0.16993091E  CI 
C.17018089E  01 
n-i7C38083E  CI 
17C56086E  01 
17C78C85E  01 
17098083E  CI 
C.17118C82F  01 
C.17138081E  CI 
0.17153079F  01 
f.l7178078E  01 
C.17193C76E  01 
0. 17218075E  CI 
C.17238073E  01 
0.17253072E  CI 
0.17273070F  CI 
f. 17298069E  01 
C.17313068E  01 
0.17338066E  01 
0.17353065E  01 
C.17373063E  01 
C.17398062E  CI 
0.1741306CE  01 
P.17433C59E  01 
0.17453057F  CI 
C.17478056F  01 
0.174C3C55F  01 
0.17518C53E  CI 
0.17538052E  01 
C.175C8050F  01 
C.17578C49E  01 
C.17^930^7E  01 
G.17fcl8046F  01 
0.17638044E  01 
C.17653043E  01 
C.17673041E  01 
C. 17693G40E  01 
O.17718039E  01 
0.17733037E  01 
C.17758036F  01 
C.17773034E  01 
0.177Q3033E  01 
C.17818031E  CI 


0.99531194F  00 
0.99581617E  00 
0.99582028E  00 
0.99582428E  00 
0.99582809E  00 
6.99583179E  00 
0.99583536E  00 
0.91583388E  00 
0.9958422°E  00 
0.9958^556E  00 
C.99584372E  00 
0.99585176E  00 
0.99535^746  00 
0.99585760E  00 
0.99586040E  00 
0.99536308F  00 
0.995«65fc4F  00 
0.99536515E  00 
0.99587053E  00 
0.99587286E  00 
0.99537512^  00 
0.99587727E  00 
0.99587935E  00 
0.99538132E  00 
0.99538323E  00 
0.995885076  00 
0.9958868CE  00 
0.99583847E  00 
0.995890C8E  00 
0.99589157E  0C 
0.99589300E  00 
0.99589437E  00 
0.99589563F  00 
0.99589694F  00 
0.99589313F  00 
C.99569Q20E  OC 
0.99590021E  00 
0.99590117E  00 
0.99590206E  00 
0.9959C29CF  OC 
0.99590367E  00 
0.99590439F  00 
0.99590504E  OG 
0.99590564E  00 
0.93590617E  00 
0.99590665E  OC 
0.99590707E  00 
0.99590749E  00 
0.99590764F  00 
0.99590814F  00 
0.99590838E  00 
0.99590«56E  00 
0.99590374E  00 
0.99590386E  00 
0.99590892F  00 
0.99590698F  00 
0.99590898E  00 
0.99590898E  00 
0.9959C898E  00 
0.99590898E  OC 
0.99590698F  00 
0.99590898E  00 
0.99590898F  00 
C.9959C898F  00 
0.99590898E  00 
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APPENDIX     II 

The  tabulated  sample  data  for  the  non-deterministic  case,    Table  II, 
was  obtained  from  a  modified  form  of  equation  (19).      Using  the 
modified  equation, 

27.  8      ,       8.  1 


0=0.  175    +    =^-2.     + 


Y(y)  X(x) 

where      X(x)    is     N(130;10) 
Y(y)     is     N(65;12) 

values  of    $    and  sensitivity  were  calculated  for  chosen  values  of  x  and 
y.     Assuming  these  to  be  measured  values,    the  data  was  then  used  to 
obtain  the  cumulative  distribution  of    0. 

Since  the  data  indicates    (f     is  more  sensitive  to  changes  in  y  than 
to  changes  in  x,    the  problem  is  worked  with  x  as  a  constant  and  y  as 
a  single  variable  as  discussed  in  section  III.     Aitkin-Hermite  inter- 
polation [Hildebrand  1956]  is  used  to  estimate  the  values  of  y  between 
data  points.      This  interpolation  scheme  utilizes  the  sensitivity  data 
and  obtains  a  more  accurate  estimate  than  would  linear  interpolation. 
The  calculation  of  the  cumulative  distribution  of    0    is  then  made 
using  equation  (35)  in  discrete  form.     As  in  the  deterministic  case, 
a      Ax  of  0.  002  was  used  for  rectangular  integration. 
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1PPENDIX  IT  -  PROGRAM  ONE 

C      CALCULATION  OF  THE  CUMULATIVE  DISTRIBUTION  OF  AN 

OUTPUT,   PHI  ,  FROM  MEASURED  VALUES  OF  THE  PARAMETERS 
X,  Yt  AND  THE  SENSITIVITY  OF   PHI   WITH  RtSPECT  TO  X 
Af^n  Y.  THE  PROBABILITY  DISTRIBUTIONS  OF  X  AND  Y  ARE 
KNOWN,  N(130:10)  AND  N(65;12)  RESPECTIVELY,  AND  THE 
DATA  POINTS  HAVE  BEEN  SELECTED  IN  DX  INTERVALS  TO 
F/CIl ITATF  THE  CALCULATIONS. 

DIMENSION  ARGl (5,5)  ,VAL(10) ,ARD  (5),XC(5),F(5,?),Z(5), 
FYT(25)  ,PXC(5)  ,X1 (200)  ,PXt 200) ,YT(25) 

C      READ  THE  MEASURED  VALUES  OF  X  AMD  THEIR  CORRESPONDING 
C      PROBABILITIES,  DENOTED  BY  XC  AND  PXC. 

READ  <  5, 102)  (X(~(I  )  ,PXC(I)  ,1  =  1,5  ) 

WPITF(6,103) 
103  FCPMAT(//,T21,  *XC«  tT47,  •PXC.  •  f /) 

WF  ITE (6,lCt)  (XC( I ) ,PXC(  I  )  ,1  =  1,5) 
106  F0PMAT(2F26.5) 

101  FCRM£T(F10.5) 

READ  THE  MEASURED  VALUES  OF  PHI  ,  OENOTED  BY  ARG1. 
THFPE  IS  A  POINT  PHI  CORRESPONDING  TO  EACH  SET  OF 
X  -  Y  POINTS. 

"t  AD  (5,101)  (  (  ARGK  I  ,  J)  ,J  =  1  ,5),  1  =  1,5) 

C      READ  THE  MEASURED  VALUES  OF  Y  AND  THE  CORRESPONDING 
SENSITIVITY  OF   PHI   WITH  RESPECT  TO  Y  FOR  EACH  OF 
THESE  VALUES,  DENOTED  BY  F(I,J). 

PEAD( 5,1021  ( (^(1  ,J)  ,J  =  1,2),I=1  ,5) 

102  F0PMATC2F10.5) 

RFAC  THE  KNOWN  CUMULATIVE  DISTRIBUTION  OF  Y,  DENOTED 
BY  YT  AMP  PYT. 

RLAD(5,1C2)  (YT(I  )  ,PYT(I  )  ,1  =  1  ,25) 
Xl( 1)=C. 5709710 
DO  20*  K=l,168 
PT=C.C 

C      X  DENOTES  THE  VALUE  OF  PHI   FOR  WHICH  THE  CUMULATIVE 
C      DISTRIBUTION  IS  BEING  CALCULATED. 

X=X1(K) 

DP  2C2  1=1  ,5 

DO  2G1  J=l,5 

Z(J)  DENOTES  THE  ARRAY  CF  ARG1  VALUES  CORRESPONDING 
TC  A  PARTICULAR  VALUE  OF  XC. 

201  7 (J)=APG1( I ,J) 
DX=1C.0 
EPS=0.CC1 
NDIM=5 

C      SUBPOUTINE  ATSM  ORDERS  THE  ARRAY  Z(J)  SUCH  THAT 
A6S(ARG(I)-X).GE.ABS(ARG(J)-X)  IF  I.GEJ.  THIS  IS 
REQUIRED  FOR  SUBROUTINE  ONE. 

CALL  ATSN(X,Z,F,5,2,ARG,VAL,NCIM) 

C      SUBROUTINE  ONE  INTERPOLATES  cUNCTION  VALUE  Y  FOP  A 
GIVEN  ARGUMENT  VALUE  X  USING  THE  TABLE  (ARG,VAL)  OF 
ARGUMENT,  FUNCTION,  AND  DERIVATIVE  VALUES  FRC"  THE 

C      SUBROUTINE  ATSM.  THIS  SUBROUTINE  USES  THE  AITKEN- 
HEFMITE  U TERPCLATICN. 
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CALL  OMP(X,ARG,VAL,Y,NDIM,EPS,IER) 


SUBROUTINE  TWO  USES  LINEAR  I  NTE  P  POL  AT  I  CN  TO  OBTAIN 
P(Y)  FPP  THc  VALUr  OF  Y  OBTAINED  IN  SUBROUTINE  ONE. 


332 

20? 

204 
107 


CAL 
PT  = 
PX( 
FOR 
CCN 
Xl( 
CCN 
WF  T 

r  c  p 

WF  I 
STO 
END 


L    TWO 

PT+(  1 

K  )=PT: 

MAK2I 

T  INI  IF 

K+1)=X1 (K)+0. 0020000 

TINUF 


'(Y,  3Y,YT,0YT) 
-FY)*PXC( I  ) 
"*DX 
1E30.8I 


TF( 


107) 
Hlf T20, 


•U( I  ) •  ,T42,' CUMULATIVE  0 I  ST R I  BUT  I  ON •  , // 


MA7(  1 

TEU  ,3021     ( XI (K)  ,PX(K)  ,K=1  ,168) 

P 


SI  BRCUTINF     1TSM  X,Z,F,  IROU.,  KCL  , APO , V AL , NC  I  w ) 

DIHFNSIPK    Z( 1 )  ,F(  1  )  ,ARG(1 ) ,7iL(  1  ) 

CACE     1^0^=1     IS    CHECKED    CUT 
IF( IPOW-l)23,21  ,1 

1  n=ntiw 

IF  N  IS  GCEATFR  THAN  IRTW,  N  IS  SET  EQUAL  TO  IRCW. 
IF(N-IfcPW)3??v2 

2  N=IPCW 


J  SUCH  THAT  Z(J)  IS  NEXT  TO  X. 


C      CASE  IRCW  .GE.  2 

C      SFARCHING  F°R  SUBSCRIPT 

3  IF(Z  (  IDpW)-Z(  1  )  )  %4,4 

4  J=IRCW 
1  =  1 

GT  Tp  6 

5  i=iptw 
J  =  i 

6  K^(J+I )/2 

IF(X-Z(K) )7,7, 8 

7  J  =  K 

GC  "T  c 
o  I=K 
-•  IF(  IAPf(  J-I  )-l  )10,10,6 

10  IF(APS(7(J)-X)-ABS(ZU)-X))1?,12,11 

11  J  =  I 

C      TABLE  SELECTION 

12  K  =  J 
JL=C 

JF=0 

DC  20  i  =  i,r 
A  P  G  ( I)  =  Z  ( K  ) 
IF( ICOL-1 )  14. 14,13 

1 3  VAL(2*I-1)=F(K) 

kk=k+ippw 

V/L( 2*1  )=F(KK) 

gc  tt  ir 

1*  VAL(  I  )=F(K  ) 
15  JJP=J+JR 

IF( JJR-IRTH)  16,18,18 
1*  JJL  =  J-JL 

IF( JJL-1  M^,  19,17 

17  IF<APS(Z(JJR+1)-X)-ARS(7(JJL-1)-X))19,19,18 

18  JL=JL+1 
K=J-JL 
GC  TO  2C 

19  JR=JP+1 
K=J+JR 

20  CONTINUE 
RFTUPN 
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CASE    IFCW=1 

21  AFG( 1 )=Z(1) 
WL(  1  )=F  ( 1  ) 

IF(ICm-2)2if22t23 

22  VAL(2)=F (2) 

23  PFTURM 
END 


C 
C 

c 


SURF  OUT INF    ONE ( X , ARC , VAL , Y , NO  I M , E  PS , I ER ) 

DIMENSION    ARG< 1 )  ,VAL(1) 

I  FP  =  2 

H2=X-£PG< 1 ) 

I F( MF I M- 1)2, 1,3 

1  Y=VAL(1)+VAL(2)*H2 

2  RETURN 

VECTOR.  APG  HAS  VORE  THAN  1  ELEMENT. 

THE  FIRST  STEP  PREPARES  VECTOR  VAL  SUCH  THAT  AITKEN 

SCHEME  CAN  RE  USED. 

3  1  =  1 

PC  5  J=2,NOIM 

H1=H2 

H2=X-£PG(J) 

Y=VAL(I ) 

VAL(  I  )=Y+VAL( I+1)*H1 

H=H1-H2 

I F  ( H  )  4  , 1 3  ,  4 

4  VAL( 1  +  1 )=Y+( VAL( I+2)-Y)*Hl/H 

5  1=1+2 

VAl  (  I  )=VAL( I  )+VAL(I+l)*H2 
END  OF  FIRST  STEP 

PREPARE  AITKEN  SCHEMF 

DELT2=0. 

IEND=I-1 

START  AITKEN-LCCP 

DC  9  I=1,IEND 

DELT1=PELT2 

Y=VAL( 1 ) 

M=(I+3)/2 

H1=APG(M) 

DC  6  J=1,I 

K=I+1-J 

L=(K+l)/2 

H=ARG(L )-Hl 

IF{H)6,14.6 

6  VAL(K)=(VAL(K) *(X-H1 )-VAL(K+l)*(X-ARG(L)))/H 
DELT2=ARS(Y-VAL(1) ) 

IF(D£LT2-EPS) 11,11,7 

7  IF( 1-5)9,8,3 

8  IF(rLLT2-DFLTl)<;t12,12 
<3  CONTINUE 

END  OF  AITKEN-LCOP 

10  Y=VAL( 1 ) 

PFTUFN 

THERt  IS  SUFFICIENT  ACCURACY  WITHIN  2*NDIM-2  ITERATION 
STFPS 

11  IFP=C 
GOTO  10 

TtST  VALUE  DELT2  STARTS  OSCILLATING 

12  IEP=1 
RETURN 

THERE  APE  TWO  IDENTICAL  ARGUMENT  VALUES  IN  VECTOR  A&G 

13  Y=VAL( 1) 

14  IER=3 
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kETUPM 
EN'D 


SIBFPUTINF    TWQ(Y,  DY,YTrPYT) 

DIMENSION    YT(25) ,PYT(25) 

IF(Y.LT.YTd)  )    GO    TC    301 

J  =  2 
110    IF(Y-YTU))     200,201,202 
202    J=J+1 

IFU.LF.25J    GO    TO    110 

PY=1.0 

GC  tc  300 
201  PY=PYT(J) 

GC  TG  30C 
20C  PY=FYT(J-1)+(PYT(J)-PYT(J-U)/(YT(J)-YT(J-1))*(Y-YT(J 

Gf  TP  3C0 
301  PY=0.0 
30C  *F  TURN 
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^TTT  CUMULATIVE 

PHI  DISTRIBUTION 

C.57C871CCE  00  0 . 62946975E-01 

0.572S7097F  00  0 . 63  1461 10F-01 

C.57437094E  00  0. 7400 1014E-01 

C.57687092E  00  0 . 79738685E-0 1 

0.57887C89E  00  0. 85873432E-01 

C.53037087E  00  0 . 92248380F-01 

C.58287C84E  00  0, 93856390E-01 

0.58487082E  00  0.1DS55303E  00 

C.58687C79E  00  0.11293674E  00 

C.58387C76E  00  0.12032640E  00 

0.59087074E  00  0.1278770CE  00 

0.59287071E  00  0.13566995E  OC 

0.594P7C69E  00  0.U363253E  OC 

0.596S7066E  00  0.15177727F  00 

C.598870fe4F  00  0.16019422^  00 

0.600870615  00  0.16903991E  00 

C.60287053E  OC  0.17780799E  00 

C.604870F6E  00  0.13684916E  OC 

G.60687053E  00  0.19600827E  00 

C.6C887051E  OC  0.20516866E  00 

C.61C87048E  00  0.21465921E  00 

0.61287045E  OC  0.22408015E  OC 

C.O1487043E  00  0.23413759F  00 

C.61687040E  00  0.24408901E  00 

C. ^18870335  OC  0.25411922E  00 

C.62C37035E  OC  0.26433289F  00 

0.62287033E  00  0.27443334E  OC 

C.6248703CE  00  0.23468227E  00 

C.62687027F  00  0.29S06922E  OC 

C.62887025E  00  0.30539531E  00 

C.63C67022E  00  0.3161C876E  00 

0.6323702CE  00  0.32671946E  00 

C.63487017E  00  0.33736047E  00 

C.63687015E  OC  0.34810561E  00 

0.63887012E  OC  0.35374993E  00 

0.640870C9F  00  0.36933064F  00 

C.64287007E  00  0.37996292c  00 

0.64467004F  00  0.390541C2E  00 

0.64687002F  00  0.40114039E  00 

C.64886999E  00  0.41177332F  00 

C.65C86997E  00  0.42229414E  00 

0.652P6994E  00  0.43286431E  00 

0.654869Q1E  OC  0.44332552F  00 

C.65636989F  00  0.45370537E  00 

C.65386986E  00  0.46397537E  00 

C.66086984E  00  0.47415489E  00 

C.66286981E  00  0.48424327E  OC 

C.66486979E  00  0.49421632F  00 

C.66686976F  00  0.50409496E  00 

0.66886973E  00  0.51387882E  00 

0.67C86971E  00  0.52356750F  00 

0.67286968F  00  0.53316504F  00 

C.67486966F  00  0.54267335E  OC 

0.67686963E  OC  0.55209064E  OC 

0.67836961E  00  0.56139952E  00 

0.68C86958E  00  0.57061529E  00 

0.68286955E  00  0.57962447E  00 

C.63486953F  00  0.58854502E  00 

0.6868695CE  OC  0.59722644E  00 

C.68886948E  OC  0.6D576999E  00 

0.69C86945E  OC  0.61422157E  00 

0.69286942E  00  0.62243212E  00 

C.6948694CE  00  0.63067019F  00 

C.69686937E  00  0.63873596E  00 

0.69886935E  00  0.64680159E  00 

0.70086932E  00  0.65472543E  OC 

C.70286930F  00  0.6625t5569E  00 

0.70486927E  00  0.67020893F  00 

C.70686924E  00  0.67763872E  00 
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C.708 
C.710 
C.712 
C.71*. 
0.716 
C.71S 
C.72C 
0.722 
0.724 
0.72  6 
0.728 
C.73C 
C.732 
0.73-^ 
0.736 
C.738 
G.740 
r.742 
0.744 
C.746 
C748 
0.750 
0.752 
C.75' 
C.  756 
0.758 
C.  76  0 
0.762 
C.764 
0.766 
C.768 
C.77C 
C.772 
C.77A 
0.776 
0.778 
C.78C 
C.732 
0.734 
C.  736 
C.788 
0.790 
r.792 
0.7^4 
C.796 
0.  793 
C.  80C 
0.  302 
0.  304 
0.  806 
C.  308 
C  .  3 1  0 
0.  512 
C.314 
0.816 
C.318 
0.  320 
C.322 
0.  324 
C.826 
C.828 
C .  8  3  C 
0.832 
0.834 
C.  336 
C.  338 
C.84C 
0.842 
0.844 
0.  846 
0.84R 
0.85  0 


86922F 
86919E 
86917E 
3691AF 
869128 
869C9E 
869C6E 
369045 
869C1C 
66  8  99E 
86896E 
36  8  9* E 
86Q91E 
868835 
36866E 
8  5 -3  93  E 
«6361E 
«63  78E 
P6«7fcP 
96373E 
8637CE 
66868E 
86365E 
86363E 
8636CE 
96858E 
86855F 
36852E 
6M5CP 
86847E 
86345E 
86342E 
r^339F 
36337E 
86«34E 
86832E 
66829F 
86827E 
86324F 
86 32  IE 
86919E 
86316E 
86314F 
86311E 
86809F 
*680^E 
363C3E 
86301F 
867°8F 
86796E 
86793F 
66791^ 
86788E 
P6785E 
86783E 
8678CE 
86773E 
F5775E 
86773R 
e6770E 
86767F 
66765E 
86762E 
36760F 
86757E 
86755E 
86752E 
86749E 
86747E 
86744E 
86742E 
86739E 


00 
00 
00 
00 
00 
00 
CO 
00 
00 
00 
00 
CO 
00 
00 
CO 

oc 

00 

oc 

00 
00 
00 

oc 

oc 

00 
CO 

oc 

00 
00 
00 
00 

oc 

00 
00 
00 
00 
00 
00 

oc 

00 
00 
CO 
00 
00 
00 
00 
00 

on 

CO 
00 
00 
OC 
CO 
00 
02 
00 
CO 
00 

cc 
oc 

00 
00 
0  0 
00 
00 
00 
00 

oc 

00 

oc 
oc 

CO 
00 


0. 

634999C4E 

OC 

0. 

6U99145F 

00 

0. 

69892645E 

OC 

0. 

70565224F 

00 

0. 

,712255845 

CO 

0. 

71330460E 

OC 

0, 

7252717CE 

00 

G. 

73166728E 

oc 

0. 

73301160F 

00 

0, 

74425936E 

00 

c. 

75045109E 

oc 

0, 

75*375795 

00 

0. 

,762222418 

oc 

0. 

75774597F 

00 

0, 

,773076^28 

00 

0, 

77835917F 

00 

0, 

78335345E 

00 

733306C0E 

oc 

0, 

,733226968 

00 

c. 

,738055538 

oc 

0. 

.30284357E 

00 

Oi 

.80759645E 

00 

0, 

81231534F 

00 

0, 

-31695616E 

00 

0, 

3215487CE 

00 

0, 

,325920105 

00 

0. 

,330205085 

00 

0, 

,83422661^ 

00 

0, 

,83805382E 

00 

c. 

,841852435 

oc 

c. 

,845386985 

00 

0, 

,848391745 

00 

0, 

,852372655 

00 

0, 

.35577905F 

oc 

0, 

,3  5  91*6858 

00 

0, 

,8625C901F 

00 

0, 

,86533793f: 

00 

0. 

.36914241F 

00 

0, 

,872419488 

oc 

0. 

.37562973E 

00 

0. 

,373737178 

oc 

0, 

,331753568 

00 

0. 

.88469028E 

00 

0, 

.33729680E 

oc 

0, 

,  83985 145E 

oc 

0, 

(/89232504E 

oc 

0, 

,83h636c1E 

00 

0, 

,83693129F 

00 

0. 

.839201C4E 

00 

0. 

•90141773E 

00 

0, 

,903618935 

00 

0. 

•9058G404E 

oc 

0. 

,90797484F 

00 

c. 

,91012955F 

oc 

c, 

,91226876E 

00 

0. 

•91439247E 

oc 

0, 

,91649936E 

00 

0, 

,918548705 

00 

c, 

,92052639E 

00 

o, 

.92240095F 

00 

0, 

,92420399E 

oc 

0. 

,925796638 

00 

0, 

,927380325 

00 

0, 

,923853755 

00 

c 

►93026102F 

00 

0. 

,931658155 

00 

0, 

.933026C8E 

00 

0. 

93436956F 

00 

0, 

,935703C2E 

00 

0, 

93702793E 

00 

0, 

.93834341E 

oc 

0. 

93964934E 

00 

64 


C.  35  2 
C.  354 

P.  856 
P.  35  8 
0.860 
P.  662 
P.  364 
C.866 
0.868 
P. 870 
P. 872 
C.874 
P.  8  76 
0.  878 
C.88C 
P.  382 
0.  334 
P.  886 
0.  338 
0.890 
C.892 
0.894 
C.  8^6 
C.  396 
C.900 
C.902 
0.904 


86736E 
66734F 
e6731E 
86729E 
86726E 
86724F 
86721E 
86718E 
86?16E 
36713E 
86711E 
86708E 
86706E 
867Q3E 
86700E 
86693E 
86695E 
86693E 
66590E 
86683E 
86 6 85 E 
86682c 
36660E 
P6677E 
36675E 
86672E 
3667CE 


00 

00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
0  0 
00 
00 
00 
CO 
00 
00 


0.94 
0.94 
0.94 
0.94 
0.94 
0.94 
0.9A 
0.94 
0.Q5 
0.95 
0.95 
0.95 
0.95 
0.95 
0.95 
0.95 
0.95 
0.95 
0.95 
0.95 
0.96 
0.96 
0.95 
0.95 
0.96 
0.96 
0.96 


094634E 
223440E 
351351E 
478309E 
603956E 
726920F 
3463C8E 
957232E 
066309E 
153637E 
250487E 
336735E 
416963E 
496714E 
575154F 
551567E 
727444F 
302903E 
377866E 
952511E 
026540E 
100211E 
173406E 
246243E 
313424E 
390605E 
461333c 


00 
OG 
00 
00 
00 
00 
00 
OC 
00 
00 
OC 
OC 
00 
00 
00 
00 
00 
00 
00 
CO 
00 
00 
00 
00 
00 
OC 
00 


65 


BIBLIOGRAPHY 


Beckmann,    P,  ,    Probability  in  Communication  Engineering,    Harcourt, 
Brace,    and  World,    1967. 

Hildebrand,    F.    B.  ,    Introduction  to  Numerical  Analysis,    McGraw- 
Hill,    1956, 

Kaplan,    W.  ,    Advanced  Calculus,    Addison- We  sley,    1957. 

Mark,    D.    G.  ,    Stember,    Jr.,    L.    H.  .    "  Variability  Analysis,  " 
Electro -Technology,    p.    37-48,    July  1965. 

Meyer,    P.    L.  .    Introductory  Probability  and  Statistical  Applications, 
Addison- Wesley,    1965. 

Nelson,    A.    C.  ,    Tommerdahl,    J.    B.  ,    "Switching  Circuit  Transient 
Performance  Prediction  Using  Empirical  Mathematical  Modeling 
Techniques",    IEEE  Transactions  on  Reliability,    p.    5-14,    March 
1965. 

Scarborough,    J.    B.  ,    Numerical  Mathematical  Analysis,    Johns 
Hopkins  Press,    1950. 

Shreider,    Y.    A.,    The  Monte  Carlo  Method,    Pergamon  Press,    1966. 

Stanton,    R.    G.  ,    Numerical  Methods  for  Science  and  Engineering, 
Prentice-Hall,    1961. 

Stiefel,    E.    L.  ,    An  Introduction  to  Numerical  Mathematics,    Academic 
Press,    1963. 

Healy,    T.    J.,    "Convolution  Revisited",    IEEE  Spectrum,    v.    6,    no.    4, 
p.    87-93,   April  1969. 

Bracewell,    R.  ,    The  Fourier  Transform  and  its  Applications, 
McGraw-Hill,    1965. 


66 


INITIAL  DISTRIBUTION  LIST 


No.    Copies 


1.  Defense  Documentation  Center  20 
Cameron  Station 

Alexandria,    Virginia  22314 

2.  Library,    Code   0212  2 
Naval  Postgraduate  School 

Monterey,    California     93940 

3.  Commander  Naval  Ships  Systems   Command  1 
Department  of  the  Navy 

Washington,    D.    C.      20360 

4.  Professor  S.    R.    Parker,    Code  52Px  1 
Department  of  Electrical  Engineering 

Naval  Postgraduate  School 
Monterey,    California     93940 

5.  Lieutenant  Alvin  L.    Franson,    USN  1 
Philadelphia  Naval  Shipyard 

Philadelphia,    Pennsylvania     19H2 

6.  Commander  1 
Naval  Electronic  Systems  Command 

Headquarters 

Washington,    D.    C.       20360 


67 


Unclassified 


Security  Classification 


DOCUMENT  CONTROL  DATA  -R&D 

(Security  classification  of  title,   body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  is   classified) 


I      ORIGINATING    A  C  T I  V  I  T  Y  (Corporal  e  author) 

Naval  Postgraduate  School 
Monterey,    California     93940 


2a.   REPORT    SECURITY    CLASSIFICATION 

Unclassified 


2b.    GROUP 


3      REPORT    TITLE 


Prediction  of  Statistical  System  Performance  from  Parameter  Distributions 


4     DESCRIPTIVE  NOTES  (Type  of  report  and.inctusive  dates) 

Master's  Thesis;  (June  1969) 


5     AUTHORIS)  (First  name,  middle  initial,  last  name) 


Alvin  LaVerne   Franson 


6      REPORT    DATE 

June   1969 


7a.     TOTAL    NO.    OF   PAGES 


67 


7b.    NO.    OF    REFS 

10 


8a.    CONTRACT    OR    GRANT    NO. 


b.    PROJEC  T    NO. 


9a.    ORIGINATOR'S    REPORT    NUMBER(S) 


9b.   OTHER   REPORT  NO(S)  (Any  other  numbers  that  may  be  assigned 
this  report) 


10.    DISTRIBUTION    STATEMENT 


Thia  aooument  has  been  approved  for  public 
r>jn4.,„-v„4.,-.  ■  . f  n.-„  ^„,irnnnt  ffi4ififfi^fvk|.sa.Te!   its  distribution  is  unlimited. 


II.    SUPPLEMENTARY    NOTES 


12.    SPONSORING   MILI  TARY    ACTIVITY 

Naval  Postgraduate  School 
Monterey,    California    93940 


13.    ABSTR  AC  T 


Techniques  for  extrapolating  statistical  data  on  the  parameters  of  individual 
components  to  the  statistical  performance  index  for  an  overall  system  are 
considered.     Two  cases  are  evaluated.     (1)    The  deterministic  case  in  which 
the   system's  performance  index  is  known  functionally  in  terms  of  the  system 
parameters.      (2)    The  non-deterministic  case  in  which  only  limited  data  on  the 
performance  index  and  its  sensitivity  with  respect  to  system  parameters  is 
known.      Computer  programs  are  developed  in  both  cases  for  combining  given 
probability  density  distributions  for  the  parameters  into  an  overall  probability 
density  distribution  for  the  system  performance  index.     Theory  and  programs 
are  developed  and  verified  with  a  specific  numerical  example. 


DD 


FORM 

1    NOV   85 

S/N    0101 -807-681 1 


1473 


(PAGE    1) 


Unclassified 


69 


Security  Classification 


A-31408 


Unclassified 


Security  Classification 


KEY     WORDS 


Cumulative  Distribution 
Deterministic 
Jacobian 

Non- deter  minis  tic 
Probability  Density- 
Sensitivity  Coefficient 


DD  ,?o":.,1473   back, 


Unclassified 


'      1  -  6  8  7  I 


70 


Security  Classification 


thesF7837 

Prediction  of  statistical  system  perform 


3  2768  000  99875  1 

DUDLEY  KNOX  LIBRARY 


